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Abstract 

In this paper we couple Synchrotron Radiation (SR) theory with a branch of phys- 
ical optics, namely laser beam optics. We show that the theory of laser beams is 
successful in characterizing radiation fields associated with any SR source. Both 
radiation beam generated by an ultra-relativistic electron in a magnetic device and 
laser beam are solutions of the wave equation based on paraxial approximation. 
It follows that they are similar in all aspects. In the space-frequency domain SR 
beams appear as laser beams whose transverse extents are large compared with 
the wavelength. In practical situations (e.g. undulator, bending magnet sources), 
radiation beams exhibit a virtual "waist" where the wavefront is often plane. Re- 
markably, the field distribution of a SR beam across the waist turns out to be strictly 
related with the inverse Fourier transform of the far-field angle distribution. Then, 
we take advantage of standard Fourier Optics techniques and apply the Fresnel 
propagation formula to characterize the SR beam. Altogether, we show that it is pos- 
sible to reconstruct the near-field distribution of the SR beam outside the magnetic 
setup from the knowledge of the far-field pattern. The general theory of SR in the 
near-zone developed in this paper is illustrated for the special cases of undulator 
radiation, edge radiation and transition undulator radiation. Using known ana- 
lytical formulas for the far-field pattern and its inverse Fourier transform we find 
analytical expressions for near-field distributions in terms of far-field distributions. 
Finally, we compare these expressions with incorrect or incomplete literature. 
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1 Introduction 



In a series of previous works [1, 2, 3] we developed a formalism which, we 
believe, is ideally suited for analysis of any Synchrotron Radiation problem. 
In particular, in [3] we took advantage of Fourier Optics ideas. Fourier Optics 
(see e.g. [4]) provides an extremely successful approach which revolution- 
ized the treatment of wave optics problems and, in particular, laser beam 
optics problems. It deals with paraxial fields, propagating along a chosen 
direction (say the longitudinal z-axis) and directions close to it. In any given 
plane transverse to the z-axis, the field is represented in the space-frequency 
domain by a complex scalar function. Alternatively, one may deal with a 
monochromatic signal oscillating as exp[—i<vt], co being the frequency and t 
being the time. Then, one can represent the slowly varying 1 field envelope 
with a complex scalar function E(r ± , z). Here z is the longitudinal position 
of the observer, and r ± is a two-dimensional vector identifying transverse 
coordinates (x,y) of the observer. E(r±,z) represents an electric field am- 
plitude polarized everywhere in the same, fixed direction. The effects of 
free-space propagation are mathematically represented by the Fresnel for- 
mula. Elementary systems, such as lenses and filters placed along the z axis 
are introduced as complex scalar-amplitude transmission functions T(r ± ), 
which yield the output field after multiplication with an input field. More 
complex systems are realized by stacking free-space sections and amplitude 
transmission functions in series. 

The use of Fourier Optics led us in a natural way to establish, in this work, 
basic foundations for the treatment of SR fields in terms of laser beam optics. 
In particular, radiation from an ultra-relativistic particle can be interpreted 
as radiation from a virtual source, which produces a laser-like beam. The 
virtual source itself is the waist of the laser-like beam, and often exhibits 
a plane wavefront. In this case it is completely specified, for any given 
polarization component, by a real-valued amplitude distribution of the field. 
As we will see, the laser-like representation of SR is intimately connected 
with the ultra-relativistic nature of the electron beam. In particular, paraxial 
approximation always applies. Then, free space basically acts as a spatial 
Fourier transformation. This means that the field in the far zone is, aside 
for a phase factor, the Fourier transform of the field at any position z down 
the beamline. It is also, aside for a phase factor, the Fourier transform of 
the virtual source. Once the field at the virtual source is known, the field 
at other longitudinal positions, both in the far and in the near zone up to 
distances to the sources comparable with the radiation wavelength, can be 
obtained with the help of the Fresnel propagation formula. This is equivalent 
to state that the near-zone field can be calculated from the knowledge of the 

1 With respect to the wavelength A = 2nc/co. 
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far-zone field, i.e. from the acceleration term of the Lienard-Wiechert field 
only. This is possible due to the paraxial approximation. On the one hand, 
as said above, the Fresnel formula can be used to recover the field up to 
longitudinal positions which are a wavelength far away from the sources. 
On the other hand, due to the paraxial approximation, the wavelength of 
interest in SR emission is very short compared with the radiation formation 
length and with the linear dimensions of the magnetic setup. In this view, it 
is incorrect to consider the far-zone field as an asymptotic expression only: 
its knowledge completely specifies, through the Fresnel integral, the near- 
zone field as well. In the case when the electron generating the field is not 
ultra-relativistic, though, the paraxial approximation cannot be applied. As a 
result, the radiation wavelength is comparable with the radiation formation 
length. It is then impossible to reconstruct the near-field distribution from 
the knowledge of the far-field pattern 2 . 

In this work we consider a SR source composed of an arbitrary number of 
magnetic devices under the standard ultra-relativistic approximation. We 
show that every such magnetic setup is equivalent to several virtual sources 
and can be modelled as a sequence of appropriately chosen plane sources 
inserted between the edges of each magnetic structure. This equivalence 
provides conceptual insight regarding SR sources and should facilitate their 
design and analysis. In fact, great simplifications in the solution for the prob- 
lem of propagation of radiation through a beamline result after adopting 
our viewpoint. Since the analysis of SR sources can be reduced to that of 
laser-like sources, it follows that any result, method of analysis or design 
and any algorithm specifically developed for laser beam optics (e.g. the code 
ZEMAX [5]) are also applicable to SR sources. 

To give an example of the applicability of our method we first consider the 
case of undulator radiation. We work within the applicability region of the 
resonance approximation. This means that we consider an undulator with 
a large number N w of periods, i.e. N w » 1, and that we are interested in fre- 
quencies near the fundamental. We find the field distribution of the virtual 
source with the help of the far-zone field distribution and we propagate to 
any distance of interest. 

Subsequently, we treat the relatively simple case of edge radiation, studying 
the emission from a setup composed by a straight section and two bending 
magnets, upstream and downstream of the straight section. Bright, long- 
wavelength edge radiation is expected to be emitted by such setup for A » 



Our arguments hold in the case near zone and formation zone coincide. For 
example, this is the case for bending magnet radiation. However, depending on the 
parameter choice, this may not be the case for complex magnetic setups. We will 
treat this subject in a more extensive way in Section 3 
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A C/ A c being the critical wavelength of SR produced by bending magnets. 
In particular, we derived an expression for the field from a straight section 
that is valid at arbitrary observation position. This result is fundamental 
importance because, due to the superposition principle, our expression may 
be used as building block for more complicated setups. We begin calculating 
an analytical expression for the edge radiation (see, among others, [6, 7, 8, 
9, 10, 11, 12, 13]) from a single electron in the far zone. Then, through 
the virtual-source technique we find novel results in the near zone. Two 
alternative description of the field propagation are given, based on a single 
virtual source located in the middle of the edge-radiation setup and based 
on two virtual sources located at the edges of the setup. 

With the help of the knowledge of the edge-radiation case, we turn to analyze 
a more complicated setup consisting of an undulator preceded and followed 
by two straight sections and two (upstream and downstream) bending mag- 
nets. We restrict ourselves to the long-wavelength region. Edge radiation 
from this kind of setup is commonly known as Transition Undulator Radi- 
ation (TUR). The first study on TUR appeared more than a decade ago [14]. 
In that work it was pointed out for the first time that, since an electron en- 
tering or leaving an undulator experiences a sudden change in longitudinal 
velocity, radiation with broadband spectral features, similar to transition 
radiation, had to be expected in the low-frequency region in addition to 
the usual undulator radiation. That paper constituted a theoretical basis for 
many other studies, developed in the years that followed. We remind here a 
few [8, 15, 16, 17, 18], dealing both with theoretical and experimental issues. 
More recently, TUR has been given consideration in the framework of large 
X-Ray Free Electron Laser (XFEL) projects. For example, a characterization 
of Coherent Transition Radiation from the LCLS undulator has been given 
in [19]. A method to obtain intense infrared/visible light pulses naturally 
synchronized to x-ray pulses from the LCLS XFEL by means of Coherent 
TUR has been proposed in [20], that should be used for pump-probe exper- 
iments in the femtosecond scale resolution. In view of these applications, 
there is a need to extend the knowledge of TUR to the near zone and to the 
coherent case 3 . 

Specification of what precedes and follows the undulator is of fundamental 
importance. As has been recognized in particular for TUR many years ago 
[8], it does make sense to discuss about the intensity distribution of TUR 
only when a detailed knowledge is given, of what precedes and follows 
the undulator in a specific setup. According to our novel approach, the two 
straight sections and the undulator will be associated to virtual sources with 
plane wavefronts. The field from the setup can then be described, in the near 



3 In this paper we limit our analysis to the single-particle case, leaving the study 
of coherent radiation for the future 
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as well as in the far field, as a superposition of laser-like beams, radiating at 
the same wavelength and separated by different phase shifts. 

Our work is organized as follows. In Section 2 we discuss our method 
in full generality. In fact, the analogy with laser beams can be extended 
to include SR from setups not discussed in this paper. In the following 
Section 3 we discuss our method in more detail and we focus on particular 
issues like the importance of the paraxial approximation, the accuracy of 
our calculations and a easily misleading relationship between velocity and 
radiation field. In Section 4 we apply our algorithm to the case of usual 
undulator radiation. Section 5 is dedicated to the study of the simplest edge 
radiation setup, composed of a straight section and two bending magnets. 
An undulator is then added to the setup in Section 6, where radiation is 
given as a superposition of laser-like beam from different virtual sources. 
Our findings are compared with results in literature in Section 7. Finally, in 
Section 8, we come to conclusions. 



2 Method for obtaining the electromagnetic field generated by an ultra- 
relativistic electron 



In the present Section 2 we discuss our method in all generality. In the next 
paragraphs 2.1 and 2.2 we show that radiation from an ultra-relativistic 
particle can be represented as a laser-like beam with a given virtual source. 
Our algorithm is based on two separately well-known expressions, the Fres- 
nel propagation formula and a far-field presentation of the electromagnetic 
field. 



2.1 Inverse problem with far-field data 



We represent the electric field in time domain E(r, t) as a time-dependent 
function of an observation point located at position r = r ± +zz = xx + yy + zz, 
where x, y and z are defined as (adimensional) unit vectors along horizontal, 
vertical and longitudinal direction. The field E(r, t) satisfies, in free space, 
the source-free wave equation: 



c 2 dt 2 



E = 



(1) 



Here c indicates the speed of light in vacuum. For monochromatic waves of 
angular frequency co, the wave amplitude has the form 
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E±(z, r ± , t) = E ± (z, r ± ) exp [-icot] + C.C. , 



(2) 



where "C.C." indicates the complex conjugate of the preceding term and 

E ± describes the variation of the wave amplitude. The vector E±_ actually 
represents the amplitude of the electric field in the space-frequency domain. 

We assume that the ultra-relativistic approximation is satisfied, that is al- 
ways the case for SR setups. In this case paraxial approximation applies [1]. 
Paraxial approximation implies a slowly varying envelope of the field with 
respect to the wavelength. It is therefore convenient to introduce the slowly 
varying envelope of the transverse field components as 

-± 

E± = E± exp [-icoz/c] , (3) 



where the propagation constant cole is related to the wavelength A by cole = 
2n/A. In paraxial approximation and in free space, the following parabolic 
equation holds for the complex envelope E of the Fourier transform of the 
electric field along a fixed polarization component: 

, 2ico d \ ~ 



The derivatives in the Laplacian operator V ± 2 are taken with respect to 
the transverse coordinates. One has to solve Eq. (4) with a given initial 
condition at z, which defines a Cauchy problem. Note that Eq. (4) is similar 
to the time-dependent Schroedinger equation. We obtain 



Ej_(z , r OJ _) — 



tco 



2nc(z - z) 



± E ± (z,r\)exp 



ico \r 01 _ - r' ± \ 
2c(z - z) 



(5) 



where the integral is performed over the transverse plane. It is important to 
note that since E is a slowly-varying function with respect to the wavelength, 
within the accuracy of the paraxial approximation one cannot resolve the 
evolution of the field on a longitudinal scale of order of the wavelength. In 
order to do so, the paraxial equation Eq. (4) should be replaced by a more 
general Helmholtz equation. 

Next to the propagation equation for the field in free space, Eq. (5), we can 
discuss a propagation equation for the spatial Fourier transform of the field, 
which can also be derived from Eq. (4) and will be useful later on. We will 
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indicate the spatial Fourier transform of Ej_(z, r' ±) with F(z, it) 



4 . 



V{z,it) = J dr'JZJ^z, f s) exp [ir' ± ■ u] . 



(6) 



Eq. (4) can be rewritten in terms of F as 

(v ± 2 + ^|_JU duF(z,u)exp[-ir ± -u\\ = 0. (7) 



Hi 2ico d \ „ , ^ 
+ — - F(z,rf) = 0. (8) 



Eq. (7) requires that 

c <9z^ 

Solution of Eq. (8) can be presented as 

ic\il\ 2 z 



F (z, u) = F (0, u) exp 



2o) 



(9) 



It should be noted that the definition of F(0, it) is a matter of initial conditions. 
In many practical cases, including the totality of the situation treated in this 
paper, F (or E) may have no direct physical meaning at z = 0. For instance, 
in all cases considered in this paper, z = is within the magnetic setup. 
However, F(0, it) can be considered as the spatial Fourier transform of the 
field produced by a virtual source. Such a source is defined by the fact 



4 For the sake of completeness we explicitly write the definitions of the two- 
dimensional Fourier transform and inverse transform of a function g(r) in agree- 
ment with the notations used in this paper. The Fourier transform and inverse 
transform pair reads: 

g(k) = J drg(r)exp[ir-k\ ; g(r) = J dk g(k) exp [-z'f • k\ , 

the integration being understood over the entire plane. If g is circular symmetric 
we can introduce the Fourier-Bessel transform and inverse transform pair: 



oo oo 

g(k) = 2n J drrg(r)J (kr); g(r) = -L J dkkg(k)J (kr) , 



r and k indicating the modulus of the vectors r and k respectively, and J being the 
zero-th order Bessel function of the first kind. 
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that, supposedly placed at z = 0, it would produce, at any distance from 
the magnetic system under study, the same field as the real source does. 
The result in Eq. (9) is very general. On the one hand, the spatial Fourier 
transform of the electric field exhibits an almost trivial behavior in z, since 
|F(z)| 2 = const. On the other hand, the behavior of the electric field itself 
is not trivial at all. These properties directly follow from the propagation 
equation for the field and its Fourier transform. 

Let us discuss the physical meaning of Eq. (9). The spatial Fourier transform 
of the field, F(z, u), maybe interpreted as a superposition of plane waves (the 
so-called angular spectrum). Once the frequency co is fixed, the wave number 
k = to I c is fixed as well, and a given value of the transverse component 

of the wave vector k ± = u corresponds to a given angle of propagation 
of a plane wave. Different propagation directions correspond to different 
distances travelled to get to a certain observation point. Therefore, they also 
correspond to different phase shifts, which depend on the position along the 
z axis (see, for example, reference [4] Section 3.7). Free space basically acts as 
a Fourier transformation. This means that the field in the far zone is, phase 
factor and proportionality factor aside, the spatial Fourier transform of the 
field at any position z. We now demonstrate this fact. This proof will directly 
result in an operative method to calculate the field at the virtual source, once 
the far field is known. We first recall that if we know the field at a given 
position (z, r' ± ) we may use Eq. (5) to calculate the field at another position 
(z ,f ±) . Let us now consider the limit z — > oo, with finite ratio r 0± /z . In 
this case, the exponential function in Eq. (5) can be expanded giving 



E±(z , r ±) ■ 



ICO 


J dr' A 


7.ncz Q 






\ico I 


x exp 


2cz \ 



\r A 2 -2r 0± -r' ± + 



I -> 12 v 

z\r ±_\ 



(10) 



Letting 6 = r OJ _/z we have 



2c 



(z + z) 



cod 



z,— 



(11) 



With the help of Eq. (9), Eq. (11) may be presented as 



Eq. (12) shows what we wanted to demonstrate: free space basically acts as 
a Fourier transformation. Also, from Eq. (12) one directly obtains 



Ex(0,r ± ) 



_ icoz r 
~2nc J 



ddexp 



ico\e\ 2 

2c 



— —T 

Ej_(d) exp 



— r ± ■ 6 

c 



(13) 



where the transverse vector r± defines a transverse position on the virtual 
source plane at z = 0. Eq. (13) allows to calculate the field at the virtual 
source once the field in the far zone is known. Identification of the position 
z = with the virtual source position is always possible, but not always 
convenient (although often it is). In general, if the virtual source is at position 
z = z s , one obtains, from Eq. (11) that Eq. (13) should be modified to 



E±(z s ,r ± ) 



_ icoz C 
~2nc J 



ddexp 



ia)\d\ : 
2c 



(z + Zs) 



Ej_(6) exp 



—r ± ■ e 

c 



(14) 



2.2 Far-field pattern calculations 



All is left to do in order to solve the propagation problem in paraxial ap- 
proximation is to calculate the field in the far zone. We indicate the electron 
velocity in units of c with jS(f'), the Lorentz factor (that will be considered 
fixed throughout this paper) with y, the (negative) charge of the electron 
with (-e) and the electron trajectory as r'(f'). Finally, the direction from the 
retarded position of the electron to the observer is fixed by the unit vector 
n, 



n = 



r'{t') 



\r ~ r'(t')\ 



(15) 



Condition n = constant defines the far-field zone. In this zone, a widely-used 
result presented in textbooks (see [21]) holds: 



-> \ la)e 
E(r 0/ a)) = exp 

cr 



ico 



-» -> 
n ■ r n 



CO 



n x (n x /3) exp 



tco 



(16) 



The transverse components of the envelope of the field in Eq. (16) can be 
written as (see [1]) 
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oo 

~ . _» , icoe C , . r ._. 
E±(Zo,^ D ,w) = — 5— dz exp[zO T 

c Z J 



'v x {z') 



x -x'(z') \ £ ^ Mz') _ y -y'{z') \ - 



where the total phase O r is 



s(z') z' 
c 







u /l z'\ [x, - x'{z ')Y + [y - y >{z>)X 

+ ii: + z!j 



2c 



(17) 



(18) 



Eq. (17) can be obtained starting directly with Maxwell's equations in the 
space-frequency domain. Textbooks (see for example [21]) usually follow a 
not so direct derivation. They start with the solution of Maxwell's equation 
in the space-time domain, the well-known Lienard-Wiechert expression, 
and they subsequently apply a Fourier transformation. Eq. (17) is automat- 
ically subject to the paraxial approximation. Here v x {z') and v y {z') are the 
horizontal and the vertical components of the transverse velocity of the elec- 
tron, while x'(z') and y'(z') specify the transverse position of the electron as 
a function of the longitudinal position. Finally, we define the curvilinear 
abscissa s(z') = vt'{z'), v being the modulus of the velocity of the electron. 

Eq. (17) can be used to characterize the far field from an electron moving 
on any trajectory as long as the ultra-relativistic approximation is satisfied. 
Then, once the far field is known, Eq. (13) (or Eq. (14)) can be used to 
calculate the field distribution at the virtual source. Finally, Eq. (5) solves 
the propagation problem at any observation position z . Note that part of 

the phase in Eq. (18) compensates with the phase in \6\ 2 in Eq. (13). If Eq. 
(17) describes a field with a spherical wavefront with center at z = 0, such 
compensation is complete. The centrum of the spherical wavefront is a 
privileged point, and the plane at z = exhibits a plane wavefront. This 
explains why the point z = is often privileged with respect to others. 

Let us conclude this Section with a short summary. We considered the 
electric field represented in the space-frequency domain. Since the system 
is ultra-relativistic, paraxial approximation applies. We saw that within the 
applicability region of the paraxial approximation, all information about 
the field, both in the near and in the far zone, is encoded in the far zone 
field alone. An explicit expression for the near field can in fact be recovered 
with the help of the paraxial Green's function for Maxwell's equation. Such 
Green's function is the same Green's function considered in Fourier Optics. 
This consideration allows to describe simple ultra-relativistic systems like, 
for instance, an electron moving through a bending magnet, an undulator or 
a straight section in terms of laser-like beams. Eachbeam has a virtual source, 
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which often exhibits a plane wavefront and is completely specified, for any 
given polarization component, by a real-valued amplitude distribution for 
the field. Note that the square modulus of such distribution can be physically 
associated with an intensity, and can be detected as the image of a properly 
placed, ideal converging lens. The Fourier transform of such distribution 
instead is linked to the electric field in the far zone. Then, once the electric 
field in the far zone is calculated, a virtual source for the laser-like beam can 
be specified. Once the field at the virtual source is known, the field at any 
position can be found by applying usual propagation equations based on 
Fourier Optics. 



3 Discussion 

In the present Section 3 we discuss the algorithm proposed in Section 2. 
First we present a thorough derivation of Eq. (17) based on paraxial Green's 
function techniques. Then we discuss the region of applicability and the 
accuracy of this derivation. Finally, we deal with some delicate question 
related to the nature of the so-called "velocity field". 

3.1 A derivation of the electric field based on paraxial Green's function 

Eq. (17) is an expression that can be used to characterize the far field from 
of an electron moving on any trajectory in paraxial approximation. For the 
sake of completeness we present, here, a derivation of Eq. (17) based on [1]. 

Accounting for electromagnetic sources, i.e. in a region of space where 
current and charge densities are present, the following equation for the 
field in the space-frequency domain holds in all generality: 

c 2 V 2 E + o) 2 E = 4nc 2 Vp-4:nia)] , (19) 

-± 

where p(r,co) and j(r,co) are the Fourier transforms of the charge density 

p(r, t) and of the current density j(r, t). Eq. (19) is the well-known Helmholtz 
equation, that has elliptic characteristic. With the help of Eq. (3), Eq. (19) can 
be written as 

c 2 exp [icoz/c] |v 2 + ^ ^- j E = 4nc 2 Vp - 4nicoJ . (20) 
A system of electromagnetic sources in the space-time can be conveniently 
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described by p(r, t) and j(r, t). In this paper we will be concerned about a 
single electron. Using the Dirac delta distribution, we can write, following 

[1] 

p(K t) = -eb{?- r'{t)) = — e —5(? ± - r\(z))S - t) (21) 

v z (z) \ V 



and 



ll t) = -et{tW- f{t)) = -^%zW± ~ r\(z))5 \^-t), 

(22) 

where r'(t) and v(t) are, respectively, the position and the velocity of the 
particle at a given time t in a fixed reference frame, and v z is the longitudinal 
velocity of the electron. 

In the space-frequency domain, Eq. (21) and Eq. (22) transform to: 

p{fl, z, co) = e —b (fl - fjzj) exp [icos(z)/v] (23) 



and 

;(rl, z, co) = e — v(z)5 (r~l - / ± (z)) exp [icos(z)/v] . (24) 

v z {z) v i 



By substitution in Eq. (20) we obtain 



x ^(z)S (rl-rl(z))-' 

Eq. (25) is still fully general and may be solved in any fixed reference system 
(x, y,z) of choice with the help of an appropriate Green's function. 

When the longitudinal velocity of the electron, v z , is close to the speed of light 
c (i.e. y\ » 1), the Fourier components of the source are almost synchronized 
with the electromagnetic wave travelling at the speed of light. In this case 
the phase co(s(z) /v -z/c) is a slow function of z compared to the wavelength. 
For example, in the particular case of motion on a straight section, one has 
s(z) = z/v z , so that co(s(z)/v - z/c) = coz/(2ylc), and if y\ » 1 such phase 
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grows slowly in z with respect to the wavelength. For a more generic motion, 
one has 



CO 



s(z 2 ) - s(zi) z 2 - Zi 



V 




2y 2 z (z)c 



CO 



(26) 



Mathematically, the phase in Eq. (26) enters in the Green's function solution 
of Eq. (25) as a factor in the integrand. As we integrate along z', the factor 
co(s(z')/v - z' '/c) leads to an oscillatory behavior of the integrand over a 
certain integration range in z' . Such range can be identified with the value 
of z 2 - Zi for which the right hand side of Eq. (26) is of order unity, and it 
is naturally defined as the radiation formation length Lf of the system at 
frequency co. It is easy to see by inspection of Eq. (26) that if v is sensibly 
smaller than c (but still of order c), i.e. v ~ c but l/yl ~ 1, then Lf ~ A. On 
the contrary, when v is very close to c, i.e. l/yl <sc 1, the right hand side of 
Eq. (26) is of order unity for Lf = z 2 - Zi » A. When the radiation formation 

length is much longer than the wavelength, E does not vary much along z 
on the scale of A, that is | d z E X/y |<sc co/c \ E XiV |. Therefore, the second order 
derivative with respect to z in the V 2 operator on the left hand side of Eq. 
(25) is negligible with respect to the first order derivative. Eq. (25) can then 
be simplified to 



where we consider transverse components of E only and we substituted v z (z) 
with c, having used the fact that l/yl <^ 1. Eq. (27) is Maxwell's equation in 
paraxial approximation. In this way we transformed Eq. (25), which is an 
elliptic partial differential equation, into Eq. (27), that is of parabolic type. 
Note that the applicability of the paraxial approximation depends on the 
ultra-relativistic assumption y 2 » 1 but not on the choice of the z axis. In 
fact, if the longitudinal z direction is chosen in such a way that yl ~ 1, 
than the formation length is very short Lf ~ A, and the radiated field is 
practically zero. As a result, Eq. (27) can always be applied, i.e. the paraxial 
approximation can always be applied, whenever y 2 » 1. Formally, this 
statement is in contradiction with the approximation l/yl « 1, explicitly 
declared to obtain Maxwell's equation in paraxial approximation, i.e. Eq. 
(27). However, suppose that one is using a detector not accurate enough 
to distinguish, on an axis where yl » 1, between solution of Maxwell's 




(27) 
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equation with paraxial approximation and solution without it. Let us fix any 
observation wavelength of interest A. Then, with the accuracy of the paraxial 
approximation, one would detect no photon flux along an axis where y\~\. 
Thus, from a practical viewpoint, the paraxial approximation can always be 
applied independently of y z . When it is not formally applicable, the error 
resulting from its application cannot be resolved by our detector, within 
the accuracy of the paraxial approximation, because we are working with 
ultra-relativistic particles, and radiation is highly collimated. In this regard, 
it should be remarked here that the status of the paraxial equation Eq. 
(27) in Synchrotron Radiation theory is different from that of the paraxial 
equation in Physical Optics. In the latter case, the paraxial approximation 
is satisfied only by small observation angles. For example, one may think 
of a setup where a thermal source is studied by an observer positioned 
at a long distance from the source and behind a limiting aperture. Only if 
a small-angle acceptance is considered the paraxial approximation can be 
applied. On the contrary, due to the ultra-relativistic nature of the electrons 
that generate radiation, Synchrotron Radiation is highly collimated. As a 
result the paraxial equation can always be applied in practice, because it 
practically returns zero field at angles where it should not be applied. 

The Green's function for Eq. (27), namely the solution corresponding to the 
unit point source, satisfies the equation 

v±2 + T^) G(z ° " z; rl ° ~ rl) = 6{rl ° ~ rl)6{z ° ~ z) ' (28) 

and, in an unbounded region, can be explicitly written as 

G(z - z'; r ±0 - f ± ) = --r—r^ jt exp — / (29) 

4n(z -z') r [ 2c(z -z') J 

assuming z — z' > 0. When z - z' < the paraxial approximation does 
not hold, and the paraxial wave equation Eq. (27) should be substituted, 
in the space-frequency domain, by the more general Helmholtz equation. 
However, the radiation formation length for z - z' < is very short with 
respect to the case z -z' > 0, i.e. there is no radiation for observer positions 
z - z' < 0. As a result, in this paper we will consider only z - z' > 0. It 
follows that the observer is located downstream of the sources. 

This leads to the solution 



oo 

E±(z 0/ r ±0/ co) = -- | dz' — — - | dr'j 
c J z - z' J 
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X 



^v ± (z')5 (r'± ~ r' ± (z'j) - V' ± 5 (r\ - r\(z')) 
x exp < z'cu 



I r ±0 - r' ± | 2 + /s(zO _ ^ 
2c(z -z') \ v c. 



(30) 



where V' ± represents the gradient operator with respect to the source point. 
The integration over transverse coordinates can be carried out leading to 
the final result: 



icoe f ^exptz'Oj] 

z -z' 



icoe C 



dz 



'v x (z') 



x -x'(z')\ Mz') y -y'{z') \ _ 
x + \ — |y 



z -z' 



z -z' 

where the total phase O r is given by 

\x -x'(z')) 2 + (y -y'(z')) 2 



(31) 



O t = co 



s(z') z' 
V c 



+ CO 



2c(z - z') 



(32) 



As we will see in Section 3.2, Eq. (31) is valid at any observation position z 
such that the paraxial approximation is valid, i.e. up to distances between 
the observer and the electromagnetic sources comparable with the radiation 
wavelength. In the far zone Eq. (31) reduces to Eq. (17). 



3.2 Exact solution ofHelmholtz equation 



Up to now we presented an algorithm based on the paraxial approximation. 
It makes sense to ask what is the range of observer positions where the 
paraxial approximation applies, and what is the accuracy of our result in 
this range. We will demonstrate that paraxial approximation holds with 
good accuracy up to positions of the observer such that its distance from the 
electromagnetic sources in the space-frequency domain, denoted with d, is 
of order of A. 

In order to investigate the accuracy of our findings we have to compare 
results from the paraxial treatment with results obtained within a more 
general framework without the help of the paraxial approximation. There- 
fore we go back to the most general Eq. (19), that can be solved by direct 
application of the Green's function for the Helmholtz equation 
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m-* "1\ 1 (■ C0 \ ~> ~1 

G(r - f) = rj — exp \ i—\ r - f 

4n \r -r'\ ^ c 



that automatically includes the Sommerfeld radiation condition 

-> as \f - r'\ — > oo . 



_ dG .cG) 



(33) 



(34) 



Integrating by parts the term in Vp we have 



■■=-!*> 



ico 



c \r - r'\ 



pn — 

c 



+ 



pn 



\r ~ r'\ 



exp 



( . co . _^ -> .) 
U—\ r -r'\ \ . 



(35) 



Use of explicit expressions for p, i.e. Eq. (23), and for i.e. Eq. (24), leads 
straightforwardly to the following expression for the field envelope: 



oo 

E{f ,co) = -— I dz'—— 
c J v z {z' 

— oo 

HP? 



n 



ic 



n 



\f - r'(z')\ M \f - r'(z')\ 2 



x exp 



r'{z')\ 



Zn ~ Z' 



(36) 



It should be emphasized that in usual Physical Optics cases one is interested 
in solving Maxwell's equations with particular boundary conditions. Typi- 
cally, boundary conditions are artificial and not consistent with Maxwell's 
equation. Thus, it is not possible to control approximations applied. On the 
contrary, in our case, Eq. (36) is an exact solution of Maxwell's equations 
with boundary conditions at infinity. In this case the scalar Green's function 
in Eq. (33) can be used without introducing any extra assumption. These 
situations are completely different. The exactness of the solution in Eq. (36) 
allows us to control the overall accuracy of the paraxial approximation 
by comparing results from Eq. (31) with results from the exact solution of 
Maxwell's equations, Eq. (36). We can present here a conservative estimate. 
Let us pose 



2n 



(37) 



When d>Lf this accuracy can be estimated to amount to ft/Lf, but quickly 
decreases as Lt » d » A remaining, at least, of order of A/d. The words 
"at least" emphasize the before-mentioned fact that this is a conservative 
estimate. 
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The most general conclusion that can be drawn from Eq. (36) is that, in- 
dependently of any approximation, the use of the space-frequency domain 
automatically presents a very natural scale, the wavelength A. As it can be 
verified by inspecting Eq. (36), the integrand term scaling as \ f - f(z')\~ 2 can 
be dropped whenever d » A. This is to be taken as a sufficient condition, 
but not necessary in general. Dropping this term is automatically enforced 
when the paraxial approximation is applied, but it is not sufficient for the 
paraxial approximation to hold. In order for the paraxial approximation to 
hold, as seen before, the system must be ultra-relativistic. 

While the present discussion about Helmholtz equation is important to 
understand the region of applicability of Eq. (31) we should stress the fact 
that the inverse field problem, which relies on far-field data only, cannot 
be solved without application of the paraxial approximation. In fact, let us 
suppose that the paraxial approximation was not applicable. Then, instead 
of Eq. (4) we should solve the homogeneous Helmholtz equation 

c 2 V 2 E + a) 2 E = . (38) 



Boundary conditions are now constituted by the knowledge of the field on 
a open surface (for example, a transverse plane) and additionally, Rayleigh- 
Sommerfeld radiation condition (separately for all polarization compo- 
nents) at infinity, i.e. 

\K-?\( d ^-i—\^Q as »oo. (39) 

\ dn co ) 

However, this is not enough to reconstruct the field at any position in space. 
In order to do so, we would need to specify the sources. In this case though, 
Eq. (35) can be directly applied. The boundary conditions specified above 
allow one to solve the direct transmission problem, but not the inverse one. 
However, as we have seen, if the paraxial approximation is applicable, the 
inverse field problem has unique and stable solution. 

Based on this discussion we can formulate a definition of several observa- 
tion zones of interest (together with their complementary zones in square 
brackets). 

• Far zone [Near zone]. Defined by d such that n = const, [for the near zone 
read: n ^ const.] throughout the integration in Eq. (36). 

• Formation zone [Non-formation zone]. Defined by d < Lf 5 [for the non- 
formation zone read: d » LA . 

5 In this region one may say that the part of the electron trajectory actively con- 
tributing to the field at the observer position is the nearest to the observer and of 
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• Radiation zone [Non-radiation zone] . Defined by d such that the electro- 
magnetic field is disentangled [for the non-radiation zone read: entangled] 
from the particle field and can be interpreted as radiation 6 . 

• 1/R-zone [l/R 2 -zone] . Defined by d such that the term in 1/R 2 = \f - 
r'{z')\~ 2 can be neglected [for the l/R 2 -zone read: cannot be neglected] in the 
integrand of Eq. (36). 

• Reconstruction zone [Non-reconstruction zone]. Defined by d such that 
a solution of the inverse problem based on far-field data can be given [for 
the non-reconstruction zone read: cannot be given]. 

In all generality, from Eq. (36) follows that the 1/R-zone always coincides 
with the reconstruction zone, because field can be reconstructed from the 
far-zone data if and only if the term containing 1/R 2 in the integrand of Eq. 
(36) can be neglected. This is clear because the far-zone data do not include 
the term in Eq. (36). Moreover if an observer is in the 1/R-zone (or in the 
reconstruction zone), it necessarily belongs to the radiation zone as well. 
However we cannot prove, in general, that if an observer is in the radiation 
zone, it necessarily belongs to the 1/R-zone (or to the reconstruction zone). 

Let us now consider a slightly less general standpoint and discuss, indepen- 
dently of the paraxial approximation, both non-ultra-relativistic systems 7 , 
such that jS 1 and y 2 - 1, and ultra-relativistic systems such that y 2 » 1 
(and obviously jS 1). 

In the case of non-ultra-relativistic systems the characteristic size of the 
system a, the formation length Lf and the radiation wavelength A are of 
the same magnitude 8 . As a result, there are only two observation zones of 
interest. The first, such that d » A, incorporates far zone, non-formation 
zone, radiation zone, 1/R-zone and reconstruction zone. The second, such 
that d < A, incorporates near zone, formation zone, non-radiation zone, 
l/R 2 -zone and non-reconstruction zone. 

In the case of ultra-relativistic systems instead, we can say that a > Lf and 
Lf » A. As a result the near zone is now defined by d < a, the formation 
zone by d < Lf, the radiation zone by d » A. As for the reconstruction zone 
and the 1/R-zone we can say that a necessary and sufficient condition to be 

length of order Lf. 

6 This means that in the radiation zone the flux of the Poynting vector through any 
transverse plane is independent of its longitudinal position. 

7 Note that there is a certain abuse of language in this definition, as it does not 
include the case jS «: 1. 

8 This is almost always the case. It is possible, however, to organize particular 
situations where contributions from many structure sum up coherently within a 
limited angular and spectral range. For instance, one may think of contributions 
from a row of bends spaced ad hoc, as in a wiggler, behaving as an array of antennas. 
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in these zones is that d » A. 



3.3 Radiation contribution from the "velocity field 



As one can see, the case of ultra-relativistic system presents an increased 
level of complexity with respect to the non-ultra-relativistic system. This 
complexity is at the origin of several misconceptions, that we here discuss. 

We presented, in the space-frequency domain, an exact solution of Maxwell's 
equation, Eq. (36) (and its paraxial form Eq.(31)), starting directly with 
Maxwell's equations in the space-frequency domain. As said before, text- 
books (see for example [21]) usually follow a not so direct derivation. They 
start with the solution of Maxwell's equation in the space-time domain, 
the well-known Lienard-Wiechert expression, and they subsequently apply 
a Fourier transformation. The Lienard-Wiechert expression for the electric 
field of a point charge (—e) reads (see, e.g. [21]): 



n-B e nx[(n- B)xB] 

E(r or t) = -e J- — - J J F i . (40) 

y 2 (l - n ■ B) 2 \r - r'\ 2 c (1 - n • B) 2 \r - r'\ 

As before, R = \r - r'(t')\ denotes the displacement vector from the retarded 
position of the charge to the point where the fields are calculated. Moreover, 

B = v/c,jl = v/c, while v and v denote the retarded velocity and acceleration 
of the electron. Finally, the observation time t is linked with the retarded 
time f by the retardation condition R = c(t - f). As is well-known, Eq. (40) 
serves as a basis for the decomposition of the electric field into a sum of two 
quantities. The first term on the right-hand side of Eq. (40) is independent of 
acceleration, while the second term linearly depends on it. For this reason, 
the first term is called "velocity field", and the second "acceleration field" 
[21]. Velocity fields differ from acceleration fields in several respects, one of 
which is the behavior in the limit for a very large distance from the electron. 
There one finds that the velocity fields decrease like R~ 2 , while acceleration 
fields only decrease as R~ l . Let us apply a Fourier transformation: 



E(r ,co) = -e 



oo 



>\2 



OO 



n-B 
y 2 (l - n- B) 2 \r - r 

dt — exp 

{l-n-B) 2 \? -r'\ 



exp 



ico t' + 



\r -r'{t'f 



ICO It H 



(41) 
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As in Eq. (40) one may formally recognize a velocity and an acceleration 
term in Eq. (41) as well. Since Eq. (41) follows directly from Eq. (40), that 
is valid in the time domain, the magnitude of the velocity and acceleration 
parts in Eq. (41), that include terms in 1/R 2 and 1/R respectively, do not 
depend on the wavelength A. On the contrary, as we have seen, the terms 
including 1 /R and 1 /R 2 in Eq. (36) do. It is instructive to take advantage of 
integration by parts to show that Eq. (36) and Eq. (40) are equivalent. With 
the help of 



-^-\K-r'(t')\ = -ft^. 
c at' 



and 



(42) 



^ = C —\-$ (43) 

dt ' \r -r'(t')\ [ P 1 PjJ 

Eq. (41) can be written as 



oo 

-± f n 

E ( r , co) = -e I dt' exp 

J \r - r\t')\ 2 

— co 

cJ dt ' [(1 -it- p)\r -r'(t') 



exp 



. L, , \r -r'(t't 

ICO It H 



(44) 



Eq. (44) may now be integrated by parts. When edge terms can be dropped 9 

9 The only assumption made going from Eq. (41) to Eq. (45) is that edge terms 
in the integration by parts can be dropped. This assumption can be justified by 
means of physical arguments in the most general situation accounting for the fact 
that the integral in dt' has to be performed over the entire history of the particle 
and that at t' - -oo and f = +oo the electron does not contribute to the field 
anymore. Let us give a concrete example for an ultra-relativistic electron. In this 
case, switchers may be obtained placing bending magnets at the beginning and 
at the end of the setup, and requiring that they deflect the electron trajectory of 
an angle much larger than the maximal observation angle of interest for radiation 
from a bending magnet. Equivalently, this means that the magnets would be longer 
than the formation length associated with the bends, i.e. Lr, ~ (cp 2 /oj) 1 ^ , where p 
is the bending radius. In this way, intuitively, the magnets act like switches: the first 
magnet switches the radiation on, the second switches it off. Then, what precedes 
the upstream bend and what follows the downstream bend does not contribute to 
the field detected at the screen position. With these caveat Eq. (45) is completely 
equivalent to Eq. (41). 
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one obtains 



oo 

E(r ,co) = I at 



$-n ic n 
K ~ r'(t')\ o) \f - r'{t')\ 2 



x 



/. L , \K-f(n\ \) 



(45) 



Recalling f = z' /v z (z') = s(z')/v and the definition in Eq. (3) one sees that 
Eq. (45) is equivalent to Eq. (36). Eq. (45) is the starting point for numerical 
codes based on direct calculation of the near field (the most well-known 
example is perhaps SRW [22]). 

The derivation of Eq. (45) is particularly instructive because shows that each 
term in Eq. (45) is due to a combination of velocity and acceleration terms 
in Eq. (41). In other words the terms in 1/R and in 1/R 2 in Eq. (45) appear 
as a combination of the terms in 1 /R (acceleration term) and 1 /R 2 (velocity 
term) in Eq. (41). As a result, one can say that there exist contributions to 
the radiation from the velocity part in Eq. (41). The presentation in Eq. (45) 
(or Eq. (36)) is more interesting with respect to that in Eq. (41) (although 
equivalent to it) because the magnitude of the l/R 2 -term in Eq. (45) can 
be directly compared with the magnitude of the 1/R-term inside the sign 
of integral, and is related, as we have seen, to the 1/R-zone and to the 
reconstruction zone, where the field can be reconstructed from far-field 
data. It must be clear that the l/R 2 -term in Eq. (45) does not coincide with 
the 1/R 2 -term in Eq. (41), that is the velocity field. If one forgets about this 
fact he would come to the paradoxical conclusion that far-field data can be 
used to reconstruct the field in the far-zone only (on the contrary, we have 
seen that when d » A we are in the reconstruction zone). Note that while 
this conclusion is, in general, a misconception, it does give incorrect results 
in the case of ultra-relativistic systems only, because in the case of non-ultra- 
relativistic systems the reconstruction zone and the far zone coincide. 

The bottom line is that only to the integral in Eq. (41) or Eq. (45) can be 
ascribed physical sense. The integrand is, in fact, an artificial construction. 
In this regard, it is interesting to note that the integration by parts giving 
Eq. (45) is not unique. An expression equivalent to Eq. (45) for the Fourier 
transform of the electric fields can be found in the very interesting, but 
perhaps little known reference [23]. After starting with Eq. (41), the authors 
of [23] organized integration by part in a different way compared with what 
has been done in Eq. (44). First they found that 
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ft x (ft x jS) 



#x[(tf-j8)Xj8] 1 
|r„ - r*|(l - n • /?) 2 |r„ - r'\ M [(l-ft-p) 



\K-r%i-fi-py 



(46) 



Note that Eq. (46) accounts for the fact that ft = (f - r'(t'))/\f* - r'(t')\ is not 
a constant in time. Using Eq. (46) in the integration by parts, authors of [23] 
obtained 



^ la)e f aA nX ( nX 
E(r ,cv) = -— dt -— — — 



ttx(rtxp) kp-rt- 2ft{n ■ p) 



x exp I ico j t' + — 



)| \r-> -r>(t>)\ 2 
r'(t')\\\ 



(47) 



Similarly as before, the edge terms have been dropped. Eq. (41), Eq. (45) and 
Eq. (47) are equivalent but include different integrands. This is no mistake, 
as different integrands can lead to the same integral. 

Another way to state this fact, in a more general way, is that Eq. (41), Eq. (45) 
and Eq. (47) can be obtained one from the other by means of the addition 
of a full derivative to the respective integrals. For example, the integration 
by parts used to transform Eq. (41) into Eq. (45) is equivalent to add to the 
integrand of Eq. (41) the full derivative dQ.il dt' of the quantity 



Qi = — 



(l-ft-pWo-r'in 



exp 



ico t' + 



r'(t't 



(48) 



while the integration by parts used to transform Eq. (41) into Eq. (47) is 
equivalent to add to the integrand of Eq. (41) the full derivative dQ, 2 /dt' of 
the quantity 

ft x {ft x /3) 

Q 2 = ^ exp 

(l-ft-(})\f -r'(t')\ 



ico t' + 



\r -r'(nV 



(49) 



There are infinite full derivatives of functions Q that can be summed to the 
integrand of the field-integral, giving infinite representations of the same 
physical quantity, the electric field. From this viewpoint, summing a full 
derivative to the integrand of the field-integral is analogous to a gauge 
transformation, that leaves invariant the observable quantity (the field), 
but transforms the electromagnetic potential in ways that may be useful in 
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particular situations. Of all possible choice of Q, the function Qi, that gives 
Eq. (45), is particularly useful because, as we have seen, the magnitude of 
the l/R 2 -term in Eq. (45) can be directly compared with the magnitude of 
the 1/R-term inside the sign of integral. However, the field obtained by 
summing full derivatives to the integrand is invariant for any choice of Q. 

Let us summarize the result of this discussion. We have seen how different 
zones scale with the other characteristic lengths of the problem, the wave- 
length A, the formation length Lf and the typical dimension of the system 
a. Complexities arise in the ultra-relativistic case, when A, a and Lf are very 
different quantities. In particular far zone, reconstruction zone and radia- 
tion zone are, in general, different concepts. Due to this fact, paradoxes may 
arise when ascribing separate meaning to velocity and acceleration part of 
the Lienard-Wiechert fields after Fourier transformation, i.e. in Eq. (40). In 
particular, we have shown that the integrand of the field-integral has no 
physical meaning. However, the representation of the field integral in Eq. 
(45) is particularly important, as it allows a better physical insight, very 
much as a particular choice of electromagnetic gauge can do in particular 
situations. 



4 Application 1. Undulator radiation as a laser-like beam 



An important exemplification of our algorithm is given here for the particu- 
lar case of undulator radiation. We assume, for simplicity, that the resonance 
condition with the fundamental harmonic is satisfied, that is 



where A w is the undulator period, while K is the undulator parameter 

K=^%, (51) 
2nm e c 2 



m e being the electron mass and H w being the maximum of the magnetic field 
produced by the undulator on the z axis. 

We also specify that position z = is taken in the middle of the undulator. 
We consider a planar undulator, so that the transverse velocity of an electron 
can be written as 

cK 

v_l(z') = sin(k w z')x. (52) 

y 
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A well-known expression for the angular distribution of the first harmonic 
field in the far-zone can be obtained from Eq. (17) or Eq. (16). Such expression 
is axis-symmetric, and can therefore be presented as a function of a single 
observation angle 6, where 

Q 2 = Q\ + Q 2 , (53) 



6 X and 9 y being angles measured from the undulator z-axis in the horizontal 
and in the vertical direction. One obtains the following distribution: 



E ± (Zo,6) = 



KcoeL u 
c 2 z y 



-Ajj exp 



.coz, 



o n2 



2c 



6' 
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o)L w 6 2 
4c 



(54) 



where the field is polarized in the horizontal direction. Here L w = A W N W is 
the undulator length and N w the number of undulator periods. Finally, Ajj 
is defined as 



Ajj = h\ 



K 2 
4 + 2K 2 , 



h (4 + 2K 2 ) ' 



(55) 



/„ being the n-th order Bessel function of the first kind. Eq.(54) describes a 
field with spherical wavefront centered in the middle of the undulator. Eq. 
(13) can now be used to calculate the field distribution at the virtual source 
yielding 



E ± (0,r 1 _) = i—A }} 

c L y 



*- 2Si (fi) 



(56) 



where Si(-) indicates the sin integral function and r± = \r±\ is the distance 
from the z axis on the virtual-source plane. Note that E ± (0,r ± ) is axis- 
symmetric. Eq. (56), that has been already presented in [3], describes a 
virtual field with a plane wavefront. Let us compare this virtual field with 
a laser-beam waist. In laser physics, the waist is located in the center of the 
optical cavity. In analogy with this, in our case the virtual source is located 
in the center of the undulator. Both in laser physics and in our situation the 
waist has a plane phase front and the transverse dimension of the waist is 
much longer than the wavelength. Note that the phase of the wavefront in 
Eq. (56) is shifted of -n/2 with respect to the spherical wavefront in the far 
zone. Such phase shift is analogous to the Guoy phase shift in laser physics. 
Finally, in laser physics, the Rayleigh range for a laser beam is presented 
in the form Zr = (co/c)w 2 , w being the radius of the beam at the location 
of the waist (i.e. at that position along z where the phase front is flat). This 
is defined, for example, by requiring that the intensity on the edge of an 
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aperture of radius w be one fourth of the intensity at the center of the radia- 
tion spot. In the undulator source case, the definition given above amounts 
to w = 0.9(cL w /co) 1/2 and z R = 0.8L W - L w . In the case of a laser beam the 
Rayleigh length is related to the resonator geometrical factor. In analogy 
with this, in the case of a undulator source the Rayleigh length is related, 
instead, to the undulator geometrical factor. 

Finally, use of Eq. (5) gives the field distribution at arbitrary observation 
position z outside the magnetic setup: 



E± (z , r ± ) = 



KcoeA 
c 2 y 



JJ 



Ei 



icor. 



-Ei 



icor. 



(57) 



where Ei(-) indicates the exponential integral function 10 . Taking square 
modulus of Eq. (56) and Eq. (57) one obtains, respectively, the intensity 
profile for the virtual source and the evolution of the intensity profile for 
undulator radiationboth in the near and in the far zone. Namely, introducing 
normalized units 

h = S^-?± , (58) 
I = (59) 



and 



z = 



(60) 



Note that the field is singular in the point z - L w /2 and r± = 0. This feature is 
related with the use of the resonance approximation. As has been seen before, the 
field distributions in the virtual source and in the far zone are linked by a Fourier 
transformation. This fact justifies the reciprocal relation linking small features in 
the near zone to large feature in the far zone (and viceversa large features in the 
near zone to small features in the far zone). Thus we may say that the singularity at 
z - L w /2 and r± - is another way of stating the well known fact that resonance 
approximation fails, in the far zone, for angles comparable with l/y z (or larger). 
From this viewpoint, singularity in Eq. (57) is not specific of our approach, but is 
intrinsically related with the use of resonance approximation. Note that, while for 
the first harmonic the far zone field does not exhibit singular behaviors at large 
angles, for the second harmonic one has a logarithmic divergence of the spectral 
flux integrated over angles of observation, due to a different behavior of the field 
distribution as a function of angles (see [24]). 
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Fig. 2. Evolution of the intensity profile for undulator radiation according to Eq. 
(62) (solid lines) and comparison with the far field asymptotic sinc 2 (0 2 /4) (dashed 
lines). The profiles according to Eq. (62) are shown as a function of angles at different 
observation distances z - 0.6, z - 0.75, z - 1.0 and z - 4.0. 

we obtain the relative intensity at the virtual source 



I(0,r ± ) = ±- 2 [n-2Si(rl)f , (61) 
and at any distance z , both in the near and in the far zone: 
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where J is defined as a normalized version of |E_J 2 in normalized units. Note 
that we use the notation 8 = r±/z in the near zone as well, because it is 
convenient, for future discussions, to present the intensity as in Eq. (62). 
Aside for scaling factor, the intensity profile in Eq. (61) can be detected by 
imaging the virtual plane with an ideal lens. 

The relative intensity at the virtual source is plotted in Fig. 1. The evolution 
of the intensity profile at different positions after the exit of the undulator 
according to Eq. (62) is plotted, instead in Fig. 2. 

To conclude, a single electron produces a laser-like radiation beam with a 

virtual source located in the center of the undulator whose transverse size 

is much larger than the radiation wavelength. Following [1], Eq. (54) and 

-> 

Eq. (56) can be generalized to the case of a particle with a given offset I and 
deflection angle if with respect to the longitudinal axis. The far-zone field, 
Eq. (54), can be generalized to 



El (z , rf,T,d) = - 



KcoeL w 

— Ajjexp 

c 2 z y 



2c 



exp 



xsinc 



coL,, 



e-fj 



4c 



(63) 



while the expression for the field at virtual source, Eq. (56), is transformed 
to: 



~ /. _» t> \ .Kcoe . 
E± (0,r/,/,r ± ) = z— A 7/ exp 

c y 



7z-2Si 



CO 



r, -I 



(64) 



The meaning of Eq. (63) and Eq. (64) is that offset and deflection of the single 
electron motion with respect to the longitudinal axis of the system result in 
offset and deflection of the waist plane. Computer codes (e.g. ZEMAX [5], 
PHASE [25]), often referred to as wavefront propagation codes have been 
developed and are now standard tools used to carry out Fourier Optics- 
based calculations. For example ZEMAX can be used for single-particle field 
calculations. Once the virtual source Eq. (64) is specified, it can be used as an 
input for computer codes. This allows propagation of the virtual wavefront 
in the presence of a complicated setup, very much like it has been used as 
an input to Eq. (5) in the particularly simple case of free-space propagation. 
One may also deal with SR from a beam of ultrarelativistic electrons. For 
example, in the case of spontaneous (incoherent) SR, radiation produced by 
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Fig. 3. Edge radiation geometry. The beam enters the system through a bending 
magnet, passes a straight section and leaves the setup through another bend. Ra- 
diation is collected at a distance z„ from the origin of the reference system, located 
in the middle of the straight section. 

an electron beam in an undulator can be treated as an incoherent collection 
of laser beams with different offsets and deflections (i.e. by summing up the 
intensities of laser-like beams with different offsets and deflections) 11 . 



5 Application 2. Edge radiation 

In this Section we apply the method described in Section 2 to the particular 
case of two dipole magnet edges. 

We restrict ourselves to the system depicted in Fig. 3 and we consider a 
single particle moving along the system. The electron enters the setup via a 
bending magnet, passes through a straight section (segment AB) and exits 
the setup via another bend. Radiation is collected at a distance z from the 
center of the reference system, located in the middle of the straight section. 
For wavelengths A » A c , A c being the critical wavelength of SR produced 
by bending magnets, the passage of the electron through the setup results 
in the emission of edge radiation (see, among the many references some 
early works and more recent developments in [6, 7, 8, 9, 10, 11, 12, 13] and 
references therein). 

In our case of study the trajectory and, therefore, the space integration in 

11 Note that if the electron beam is distributed coherently, radiation can be described 
as a coherent collection of laser beams. 
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Eq. (17) can be split in three parts: the two bends, that will be indicated with 
bi and b 2 , and the straight sections AB. One may write 

^± ^± ^± ^± 

E±(z , r ±0 , co) = E bl (z , r ±0 , co) + E AB (z , r ±0 , co) + E b2 (z , r ±0 , co) , (65) 
with obvious meaning of notation. 

We will denote the length of the segment AB with L. This means that points 
A and B are located at longitudinal coordinates z A = -L/2 and z B = L/2. 

First, with the help of Eq. (17) we will derive an expression for the field in 
the far zone. The intensity distribution in this case, result to be in agreement 
with expression in references cited above. Then, we will calculate the field 
distribution at the virtual source with the help of Eq. (13). Finally, Eq. (5) 
will allow us to find an expression for the field both in the near and in the 
far zone. 

In the following we will only deal with a contribution of the electric field, 
i.e. that from the straight section. 

According to the superposition principle, one should sum the contribution 
due to the straight section to that from the bends. In some cases one can 
ignore the presence of the bending magnets with good accuracy, and treat 
them as if they had zero length. In other cases one cannot do that. When 
one can neglect the field contribution due to the switchers, one can work 
in what can be called the "zero-length switchers approximation". While we 
direct the reader to Section 5.3 for details, we can mention here that for long 
enough straight sections L » y 2 A, a condition for the zero-length switchers 
approximation to apply is A » A c , A c being the critical wavelength for SR 
from bends. Intuitively, magnets act like switchers: the first magnet switches 
the radiation harmonic on, the second switches it off. In the case depicted in 
Fig. 3 the switchers are bending magnets, but other setups can be considered 
where they have different physical realizations. For example, one may think 
of a setup where the phenomenon at study is the emission of bremsstrahlung 
in a collision between a ultra-relativistic electron and a nucleus. In this case, 
the switch-off process is taken care of by the collision itself. The nucleus 
plays the role of the switcher and the impact parameter, i.e. the minimal 
distance of the nucleus to the electron, characterizes the switcher itself. 
Note that there is no principle nor practical limitation to the length of the 
switcher. In the bremsstrahlung case we can assume that such length is very 
short, even with respect to the radiation wavelength. In the case of bends it 
can be much longer than all characteristic scales of interest, depending on 
the magnetic field strength, and all kind of intermediate situations can be 
realized. As one can see, the case of a straight line preceded and followed 
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by switchers has a number of physical realizations. The only feature that 
different realizations must have in common by definition of switcher is 
that the switching process depends exponentially on the distance from the 
beginning of the process. Then, a characteristic length d s is associated to 
any switcher. Since electrodynamics is a linear theory in the case when 
the contribution from the finite straight-line cannot be considered a good 
approximation of the total field, it can be considered as a building block 
for a more complicated setup. An example of a more sophisticated setup is 
studied in the next Section 6, where the problem of Transition Undulator 
Radiation is addressed. 

In general, one can say that the problem of studying the field from a fi- 
nite straight-line is of fundamental importance, independently of the type 
of switchers considered. This justifies our attention to the straight-section 
contribution to the field. 



5.1 Far field pattern of edge radiation 

5.1.1 Field contribution calculated along the straight section 
Accounting for the geometry in Fig. 3 we have 

s(z') = z' for z A < z' < z B (66) 



With the help of Eq. (17) we write the contribution from the straight line AB 
as 

L/2 

Eab = ~y~ f dz'exp[i<I> AB ](9j+9 v $) (67) 
c z J 

-L/2 



where O^g in Eq. (67) is given by 



2c 



* + |(£ +fl 2 +fl 5) 



(68) 



6 X = x /z and 6 y = y /z being the observation angles in the horizontal and 
vertical direction. From Eq. (67) one obtains 



— icoeL 



icod 2 z 
2c 



6 sine 



4c 



r 



(69) 
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where, as before 

e 2 = 6 2 x + e 2 y 



(70) 



Eq. (69) is an exact expression for the electric field from the straight section 
AB. Note that Eq. (69) describes a spherical wave. Moreover, it explicitly 
depends on L (this last remark will be useful later on). 

The formation length Lf s for the straight section AB can be written as 

Lfe ~ min[y 2 A,L] . (71) 



Depending on the wavelength of interest then, L fs ~ y 2 A or L fs ~ L. In both 
cases, with the help of Eq. (68) we can give on a purely mathematical basis 
an upper limit to the value of the observation angle of interest related to the 
straight line: 

el, y < 4 • (72) 



Note that, if L fs ~ y 2 A, the maximal angle of interest is independent of the 
frequency. 

Finally it should be remarked that the the far-zone asymptotic in Eq. (67) 
is valid at observation positions z » L. This is a necessary and sufficient 
condition for the vector n pointing from the retarded position of the source 
to the observer, to be considered constant. This result is independent of the 
formation length given by Eq. (71). Therefore, when L < y 2 A we can say 
that an observer is the far zone if and only if it is located many formation 
lengths away from the origin. This is no more correct when L » y 2 A. In this 
case the observer can be located at a distance z » y 2 A, i.e. many formation 
lengths away from the origin of the reference system, but still at z ~ L, i.e. 
in the near zone. As we see here, the formation length Lf is often, but not 
always related with the definition of the far (or near) zone. This is the case for 
bending magnet radiation, but not for edge radiation. As has been remarked 
before in Section 3, the far (or near) is related with the characteristic size of 
the system a (in our case a = L). In its turn Lf < a, which includes, as in the 
edge radiation case when y 2 A <sc L, the situation Lf <<c a. 

5.1.2 Energy spectrum of radiation 

The radiation energy density as a function of angles and frequencies can be 
written as 
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dD. being the differential of the solid angle Q. Here we have used Parseval 
theorem and the fact that the irradiance (i.e. the energy per unit time per 
unit area, averaged over a period of the carrier frequency) is given by 
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-> 
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/ 2n 





(74) 



where brackets < ... > denote averaging over a cycle of oscillation of the 
carrier wave. Substituting Eq. (69) in Eq. (73) it follows that 



dW _ co 2 e 2 E 2 
dcodD, ~ 47i 2 c 3 



6 2 sine 2 



4c 



(75) 



It is convenient to introduce normalized quantities 




2 coE 

9 = ~r ' 

y 2 c 



(76) 
(77) 



With Eq. (75) in mind and using normalized units, we may write the direc- 
tivity diagram I of the radiation as 



I = const. X 6 2 sine 2 - (d 2 + 



(78) 



The directivity diagram in Eq. (78) is plotted in Fig. 4 for several values 
of (p as a function of the normalized angle 6. The natural angular unit is 
evidently (2ttL/A)- 1/2 . 

There are two asymptotic cases for the problem parameter ^ « 1 and 
(£> » 1. When <^> » 1 the oscillating structures are fine with respect to the 
envelope of the directivity diagram. The intensity peaks where the envelope 

of the directivity diagram peaks, i.e. at 9 ~ -yj$>. This feature for the far-field 

emission is shown in Fig. 4 for the case § = 50 and are typical for » 1. 
When the length of the straight section becomes smaller and (£> <sc 1, we 
reach the other asymptotic limit. In this case, the period of the oscillations 

due to the sin(-) function becomes much larger than J§, as it can also be 
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Fig. 4. Directivity diagram (solid lines) of the radiation from the setup in Fig. 3 for 
different values of <p. Dotted lines show the envelope of the directivity diagram. 

seen in Fig. 4 for the case (p = 0.1. Then, the maximum intensity does not 
coincide anymore with the peak of the envelope of the directivity diagram, 
but it is found at = 2.2. 

The behavior of the far-field emission described here is well-known in liter- 
ature, and has been first reported long ago in [8]. We take this as the starting 
point for further investigations based on Fourier Optics. 



5.2 Method of virtual sources 

5.2.1 Edge radiation as a field from a single virtual source 

Eq. (13) and Eq. (69) allow one to characterize the virtual source through 



E(0,rl) = 



to 



2 eL 



2nc 3 



dO sine 



4c 



exp 



10)^ * 

— r ± ■ 

c 



(79) 
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The Fourier transform in Eq. (79) is difficult to calculate analytically in full 
generality However, simple analytical results can be found in the asymptotic 
case for <ft «: 1, i.e. for 2nL/(y 2 A) «: 1. In this limit, the right hand side of 
Eq. (79) can be calculated with the help of polar coordinates. An analytic 
expression for the field amplitude at the virtual source can then be found 
and reads: 



E(0, r\) = |pxsinc g |r~l| 2 ) . (80) 

It is useful to remark, for future use, that similarly to the far-field emission 
Eq. (69), also the non-normalized version of the field in Eq. (80) explicitly 
depend on L. This is true in general, for any value of (p. After definition of 
the normalized transverse position 




the intensity distribution of the virtual source is given by 



I(f ± ) = const, x r 2 sine 2 , (82) 

it can be detected (aside for scaling factors) by imaging the virtual plane 
with an ideal lens, and is plotted in Fig. 5. 

Note that Eq. (80) describes a virtual source characterized by a plane wave- 
front. Application of the Fresnel propagation formula, Eq. (5) to Eq. (80) 
allows one to reconstruct the field both in the near and in the far region. We 
obtain the following result: 



E(z , r ± ) = T exp 

c n ± 



con 



x 



2cz„ 



icon 



icon 



exp | 



2cz (l+2z /L) / 



exp | 



2cz (-l + 2z /L). 



(83) 



Eq. (83) solves the field propagation problem for both the near and the far 
field in the limit for cp <<c 1 12 . 



12 Eq. (83) is singular as r±_ - and z = L/2. Moreover, the spectral photon 
flux, integrated in angles, is logarithmically divergent. In the undulator case, as 
we have seen in footnote 10, singular features have been interpreted as limits 
of applicability of the resonance approximation. Similarly we can state here that 
the paraxial approximation is valid, in the far zone, for angles of observation much 
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Fig. 5. Intensity distribution at the virtual source, I, as a function of \f±\ (upper plot) 
and 3D view as a function of x and y. 
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The intensity profile associated with Eq. (83) is given by 



e 2 



exp| 



id 2 z 



2(l+22 ), 



exp| 



Z0 2 Z o 



2(-l+2z ), 



(84) 



where = f_ L /z and z = z /L. 

In the limit for z — > oo Eq. (83) and Eq. (84) respectively transform to 



— icoeL 
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(85) 



and 



I = const, x 9 sine 



9 2 



(86) 



corresponding to Eq. (69) and Eq. (78) in the asymptotic case for for <ft <<c 1. 

Note that when § «; 1 we have only two asymptotic regions with respect 
to z : the far zone for z » 1 and the near zone for z <sc 1. Of course, it 
should be stressed that in the case z «: 1 we still hold the assumption 
that the approximation of zero-length switchers is satisfied. It is interesting 
to study the evolution of the intensity profile for edge radiation along the 
longitudinal axis. This gives an idea of how good the far field approximation 
(z » 1) is. A comparison between intensity profiles at different observation 
points is plotted in Fig. 6. 

The case (J « 1 studied until now corresponds to a short straight section 
L «: y 2 A. When this is condition is not satisfied, we find that the integral on 
the right hand side of Eq. (79) is difficult to calculate analytically. However, 
one can use numerical techniques. With the help of polar coordinates, the 
right hand side of Eq. (79) can be transformed in a one-dimensional integral, 
namely 



smaller than unity. The logarithmic singularity of the flux and the singular behavior 
at f ± = and z - L/2 correspond to a logarithmic divergence of the flux in the far 
zone at large angles. These divergences, in the near as well as in the far zone are 
outside the region of applicability of our approximation, because we discuss the 
far-zone field within angles much smaller than unity as well as near-zone field for 
transverse displacements r±_ » A, where we have no singularity at all. 
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Fig. 6. Evolution of the intensity profile for edge radiation in the limit for § <?C 1. 
These profiles are shown as a function of angles at different observation distances 
z - 0.6, z - 1.0, z - 2.0 and z - 5.0 (solid lines). The dashed line always refers 
to the far-zone intensity. 



E(0,A) = - 
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CO 
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e 2 
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4c 



(87) 



We calculated the intensity distribution associated with the virtual source 
for values § = 0.1, § = 1, <£> = 10 and ^ = 50, corresponding to directivity 
diagrams in the far zone in Fig. 4. We plot these distributions in Fig. 7. 
In particular, it is instructive to make a separate, enlarged plot of the case 
<ft = 50, that is in the asymptotic case for <ft » 1. This is given in Fig. 8. 
Fine structures are now evident, and are consistent with the presence of fine 
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Fig. 7. Intensity profiles of the virtual source for the setup in Fig. 3. These profiles 
are shown for (p - 0.1, (p - 1.0, (p - 10.0 and <ft — 50.0 (solid lines). Solid curves 
are calculated with the help of Eq. (87). The dotted lines show comparison with the 
asymptotic limit for § <SC 1, shown in Fig. 5 and calculated using Eq. (82). 



structures in Fig. 4 for the far zone. In fact, as we have seen, the field in the 
far region and the virtual source are linked, basically, by a Fourier transform. 
In principle, this allows to qualitatively describe the situation at the virtual 
source through the reciprocal relation. However, as we have said before, 
when (p » 1, in the far zone we have fine structures with variable width 56 
(see also Fig. 4). As a result, in the limiting case § » 1, use of the reciprocal 
relation to describe the properties of the virtual source is problematic. For 
example, a typical width 56 ~ 1 should correspond to typical dimension 
of the virtual source of order unity, that is in obvious disagreement with 
Fig. 8. Nonetheless, we managed to specified the field at the virtual source 
by means of numerical techniques, even in the case § » 1 (see Fig. 8). 
Once the field at the virtual source is specified for any value of ^, Fourier 
Optics can be used to propagate it. In free space, the Fresnel formula must 
be used. However, we prefer to proceed in another way There is, in fact, 
an alternative way to obtain the solution to the field propagation problem 
valid for any value of <ft and capable of giving a better physical insight for 
large values of 
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Fig. 8. Intensity distribution at the virtual source for the setup in Fig. 3 for <p = 50 
(enlargement of the bottom right graph in Fig. 7). 

5.2.2 Edge radiation as a superposition of the field from two virtual sources 
Let us begin considering the far field in Eq. (69). This can also be written as 



E(z ,6) = E 1 (z ,e) + E 2 (z ,e) (88) 
where 



E h2 (z ,9) = ± 



led 
cz (d 2 + 1/y 2 ) 



exp 



icoL 
~4cy 2 



exp 



icohd 2 
2c 



(89) 



The two terms E\ and E 2 represent two spherical waves respectively cen- 
tered at z = L/2 and z = -L/2, that is at the edges of the straight section. 
Analysis of Eq. (89) shows that both contributions to the total field are 
peaked at an angle of order 1/y. While, as has been seen before, the total 
field in dimensional units explicitly depends on the straights section length 
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E, the two expressions E\ and E 2 exhibit dependence on E through phase 
factors only. This fact will have interesting consequences, as we will discuss 

later. The two spherical waves represented by Ei and E 2 may be thought 
as originating from two separate virtual sources located at the edges of the 
straight section. One may then describe the system with the help of two 
separate virtual sources, and interpret the field at any distance as the super- 
position of the contributions from two edges. This viewpoint is completely 
equivalent to that considered before involving a single virtual source in the 
straight line center. We are presenting here a different description of the 
same phenomenon. As we have seen before, we could not specify, analyti- 
cally, the single virtual source in the center of the straight line. In contrast to 
this it is possible to specify the two virtual sources at the edges of the setup. 
In order to so so we take advantage of a slightly modified version of Eq. (13) 
that accounts for an arbitrary position of the source 25(1,2): 




(90) 

-> -> 

Separately substituting Ei and E 2 into Eq. (90), and with the help of polar 
coordinates, we find the following expressions for the field at the virtual 
source positions z sl = E/2 and z s2 = -E/2 13 : 



vlJ- - \ - 2iea} 



13 It should be noted here that the virtual sources are singular in the point r± - 0, 
due to the behavior of the modified Bessel function K\ . In footnote 10 we associated 
the singularity of the undulator field with the resonance approximation. Similarly 
we can state here that paraxial approximation is valid, in the far zone, for angles 
of observation much smaller than unity. As a result, features near the downstream 
edge of the straight section cannot be resolved in paraxial approximation because 
they depend on far field data at large angles. Analysis of Eq. (91) shows that the 
intensity associated to each source exhibits a weak, logarithmic singularity of the 
flux in the near zone around r± — 0. This corresponds to a logarithmic divergence 
of the flux in the far zone at large angles. These divergences, in the near as well as in 
the far zone are outside the region of applicability of our approximation, because 
we discuss the far-zone field within angles much smaller than unity as well as 
near-zone field for transverse displacements r± » A, where we have no singularity 
at all. 
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where Ki(-) is the modified Bessel function of the first order. Analysis of Eq. 
(91) shows a typical scale related to the sources dimension of order cy/co in 

dimensional units, corresponding to 1/ -yj$> in normalized units. This is in 
agreement with the fact that both source contributions to the far field are 

peaked at an angle of order Note that this remark is only qualitative. 
The peak angle in the far zone does not correspond univocally to the width 

of the field distribution, nor the typical scale 1/ of the virtual sources 
can be univocally associated to the width of the sources. Also note that 
since the far zone field in dimensional units exhibit dependence on L only 
through phase factors only has its counterpart in the fact that the field at the 
virtual sources, written in dimensional units, exhibit dependence on L only 
through phase factors as well. 

Application of the Fresnel formula allows to calculate the field at any dis- 
tance z in free space. Of course, Eq. (91) can also be used as input to any 
Fourier code to calculate the field evolution in the presence of whatever 
optical beamline. However, we restrict ourselves to the free-space case. 
Taking advantage, once more, of polar coordinates and using the definition 

6 = f ± /z we obtain the following scaling law for the intensity in normalized 
units: 
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(92) 



In the limit for z » 1 Eq. (92) gives back the square modulus of Eq. (89) 
in normalized units. Similarly, in the limit for ^ « 1, and using the fact 

that -^>fj_) - l/(fj_ yffy) one recovers the square modulus of Eq. (83). In 
general, the integrals in Eq. (92) cannot be calculated analytically, but they 
can be integrated numerically. 
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$ = 50.0 



Fig. 9. Intensity profile for edge radiation for (p = 50. These profiles are shown as a 
function of angles at different observation distances z = 0.52, z - 0.6, z - 1.5 and 
z = 100.0. 

First we checked that we are able to recover, posing z = 0, the intensity 
profile for the single virtual source, already calculated numerically and 
shown in Fig. 7 and Fig. 8. Then we propagated the field at non-virtual 
positions, for z > 1/2. In Fig. 9 we plotted, in particular, results for the 
propagation in case (p = 50. Radiation profiles are shown as a function of 
angles 9 at different observation distances z = 0.52, z = 0.6, z = 1.5 and 
z = 100.0. 

From a technical viewpoint, it is easier to deal with two sources than with 
one, because the expression for the virtual two sources is analytical, whereas 
that for a single one is not. Moreover, as said before, the two-sources picture 
gives new physical insight for the asymptotic limit <p" » 1. 

Qualitatively, we can deal with two limiting cases of the theory, the first for 
cp* «c 1 and the second for § » 1. 

Let us first discuss the case (p" <sc 1. The field at any observation distance is 
given by Eq. (83). There are only two observation zones of interest. 

• Far zone. In the limit for z » 1 one has the far field Eq. (85). 
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• Near zone. When z < 1 instead, one has the near field Eq. (83). 

As it can be seen from Eq. (83), the total field is given, both in the near 
and in the far zone, by the interference of the two virtual sources located at 
the straight section edges. Eq. (92) shows that the transverse dimension of 
these virtual sources is given by yX in dimensional units. This is the typical 
scale in after which the integrands in dr' ± in Eq. (92) are suppressed by 
the function K\. Thus, the sources at the edges of the straight section have 
a dimension that is independent of L. In the center of the setup instead, the 
virtual source has a dimension of order VAL as it can be seen Eq. (82). When 
(£> «: 1 the source in the center of the setup is much smaller than those at the 
edge. This looks paradoxical. The explanation is that the two contributions 
due to edge sources interfere in the center of the setup. In particular, when 
« 1 they nearly compensate, as they have opposite sign. As a result of 
this interference, the single virtual source in the center of the setup (and its 
far-zone counterpart) has a dimension dependent on L (in non-normalized 
units) while at the edges (and in their far-zone counterpart) the dependence 
on L is limited to phase factors only. Due to the fact that edges contributions 
nearly compensate for « 1 one may say that the single-source picture is 
particularly natural in the case (p <sc 1. 

Let us now discuss the case <£> » 1. In this situation the two-sources picture 
becomes more natural. Let us define with rf 1/2 = z + L/2 the distances of 
the observer from the edges. As seen before, the transverse dimension of 
the sources at the edges of the straight section is r' ± ~ yX. Moreover, we see 
from Eq. (67) and Eq. (68) that when <ft » 1 the formation length is Lf = y 2 X, 
much shorter than the system dimension L. As a result, one can recognize 
four regions of interest, that are more naturally discussed in the two-source 
picture. 

• Two-edge radiation. Far zone. When d i/2 » L we are summing far field 
contributions from the two edge sources. This case is well represented 
in Fig. 9 for z = 100, where interference effects between the two edges 
contribution are well visible. 

• Two-edge radiation. Near zone. When d 1/2 ~ L the observer is located 
far away with respect to the formation length of the sources L',. Both 
contributions from the sources are important, but that from the nearest 
source begins to become the main one, as &\ and d 2 become sensibly 
different. This case is well represented in Fig. 9 for z = 1.5. 

• Single-edge radiation. Far zone. When y 2 X <sc d\ <sc L the contribution 
due to the near edge is dominant, while the far edge contribution is negli- 
gible. Such tendency is clearly depicted in Fig. 9 for z = 0.6. Interference 
tends to disappear as the near edge becomes the dominant one, while the 
intensity distribution tends to approximate 
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where E 2 is the single edge far field limit in Eq. (89). 
• Single-edge radiation. Near zone. When < d\ < y 1 K we have the near- 
field contribution from a single edge. As d\ becomes smaller and smaller 
the intensity distribution tends to reproduce the singular behavior from 
a single virtual source, i.e. the square modulus of Eq. (91): 

I {^.'^ ± ) = constx 

This tendency can be seen in Fig. 9 for z = 0.52. 

It should be noted that the straight section contribution has been calculated 
within the applicability region of the paraxial approximation, but indepen- 
dently of all assumptions on the switchers. If switchers or other setup parts 
are present, their contribution must be separately calculated and added to 
the contribution due to the straight section. For instance, if an observer is 
placed after a switcher located downstream of the straight section, the field 
due to the straight section has physical meaning (i.e. the theoretical inten- 
sity can be compared with experiment) only if the zero-length switchers 
approximation can be applied. Otherwise it is just a mathematically conve- 
nient term, a partial result to be summed with the switcher contribution. 
If, instead, the observer is placed on a transverse plane placed before the 
switcher, at the downstream edge of the straight section, the switcher (or 
any other part of the setup following the straight section) will not influ- 
ence the field with ultra-relativistic accuracy. Eq. (91) (or Eq. (94)) can then 
be compared with experimental results. As has been discussed in footnote 
13, the singularity in the field at r± = is fundamentally related with our 
ignorance about the structure of the electron. 

It should be noted that our Eq. (91) is identical to the frequency-domain 
expression for the transverse component of the field originating from an 
ultra-relativistic electron moving with constant velocity (see [21]). In fact, 
our treatment for the straight-section contribution in the case of single-edge 
radiation in the near zone gives, within an accuracy of 1/y 2 , the same result 
found in [21] for bremsstrahlung in a collision between an ultra-relativistic 
electron and an atomic nucleus by means of the method of virtual quanta. 
In the bremsstrahlung case the electron harmonic is switched off within a 
length d s <sc A even in the model for an electron on a uniform motion, be- 
cause non-negligible contributions to the field are generated by the part of 
the trajectory before the nucleus only, while the trajectory after the nucleus 
gives negligible contribution within an accuracy 1/y 2 . Exact calculation, 
retaining all parts of the Lienard-Wiechert field show that the field from 



= const x 



4e 2 



y 2 /t 2 c 2 




(94) 
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the relativistic particle is equivalent to two virtual pulses of radiation. One 
moves longitudinally, i.e. in the same direction as the electron. A second 
moves transversely, i.e. perpendicularly with respect to the electron direc- 
tion. As already said, the component of the frequency spectrum due to the 
longitudinal pulse is given by our Eq. (91). At any observation position, the 
component due to the transverse pulse can be neglected with an accuracy 
1/y 2 . In other words, the weak transverse pulse is a correction to the main 
longitudinal pulse. It may be interpreted as an evanescent wave decaying in 
the longitudinal direction, since it propagates transversely, i.e. at large angle 
n 12 with respect to the z axis. Dropping this term gives back our result, that 
was found in the paraxial approximation. 

Our previous result in Eq. (91) constitutes a novel finding in theory of edge 
radiation dealing with the situation of a particle in straight motion. The case 
of a particle in straight motion is of fundamental importance, because of 
two reasons. First, our result can be applied to any setup where the zero- 
length switchers approximation applied. Second, straight sections are basic 
building blocks for any magnetic setup, and our expression for the near 
field from a straight section can be used, as part of the total field, also in the 
case when the zero-length switchers approximation fails. 



5.3 Supplementary remarks on the zero-length switcher approximation 

For the sake of completeness it remains to qualitatively discuss the condi- 
tions when the field due to switchers is negligible, because in this case one 
can work in the "zero-length switchers approximation". We consider the 
case when switchers are present in the form of bending magnets. In par- 
ticular, here we report results from an analysis of the problem in the limit 
L » y 2 A. In this case the formation length associated with the straight sec- 
tion is given by y 2 A, while the formation length associated with the bending 
magnet is given by (Ap 2 ) 1 ^ 3 , p being the radius of the bend. The ratio be- 
tween the latter and the former is indicated with e 2 = (A c /A) 2 ^ 3 <c 1, where 
A c = 4np/(3y 3 ) is the critical wavelength for SR from the bend and we as- 
sumed A » A c , as we are interested in edge radiation. The product between 
the small quantity e and the formation length for the straight section is 
indicated with d c = ey 2 A, and turns out to constitute an extra characteristic- 
length for our system. Let us indicate with d the distance between a trans- 
verse plane of observation and the downstream edge of the straight section. 
At very short distances, for d <<c d c , the straight section contribution domi- 
nates the bending magnet one for transverse displacements from the z-axis 
r± of order r± > eyA. Within the area r± < eyA the two contributions are 
comparable. Note that in this area the contribution from the straight section 
is much larger than in the area for eyA < r ± < yA due to singular behavior 
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Screen 




Fig. 10. Edge radiation geometry. The beam enters the system through a bending 
magnet, passes a straight section, an undulator and another straight section before 
leaving the setup through another bend. Radiation is collected at a distance z from 
the origin of the reference system, located in the middle of the undulator. 

of the field contribution from the straight section at r± — > 0. When d ~ d c 
the zero-length switcher approximation cannot be used, and contributions 
from the bend should be calculated explicitly. When y 2 A > d » d c we have 
room for direct application of our theory for the straight section within the 
range < r± <sc (A 2 p) 1/3 , because the characteristic size of the radiation is 
e(A 2 p) 1/3 . Finally, for d » y 1 K, we have room for direct application of our 
theory for the straight section within the angular range < 6 <sc (A/p) 1/3 , 
because the characteristic angle of the radiation is e(/t/p) 1/3 . 



6 A more sophisticated application. Transition undulator radiation 



In this Section we apply the method described in Section 2 to the more 
complicated case of an undulator setup. 

Instead of the setup in Fig. 3, we now consider the system depicted in 
Fig. 10 and we consider a single particle moving along the system. The 
electron enters the setup via a bending magnet, passes through a straight 
section (segment AB), an undulator (segment BC), and another straight 
section (segment CD). Finally, it exit the setup via another bend. Radiation 
is collected at a distance z from the center of the reference system, located in 
the middle of the undulator. The passage of the electron through the setup 
results in collimated emission of radiation in the low photon energy range, 
a mechanism analogous to transition radiation. This kind of radiation is 
known in literature as TUR [14, 8, 15, 16, 17, 18], even though in this paper 
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we prefer to denote it as edge radiation from undulator setup. 

In our case of study the trajectory and, therefore, the space integration in 
Eq. (17) can be split in five parts: the two bends, that will be indicated with 
b\ and b 2 , the two straight sections AB and CD and the undulator BC. One 
may write 



E±(z , rj_ 0/ a)) - E bl (z 

or f±o/ 

oo) + E AB (z oo) + E BC (z 

or T±o/ 

00) 

-» -> 

+E CD (z , f ±0 , oj) + E bl {z , r ±0/ oj) , (95) 
with obvious meaning of notation. 

We will denote the length of the segment AD with L tot , while we will indicate 
the length of the straight section AB with L\, the length of the straight section 
CD with L 2 and the length of the undulator with L w . It follows 

L tot = U+ Lw + L 2 • (96) 



This means that point A is located at longitudinal coordinate Za = —L\—L W I2, 
while B, C and D are located respectively at z B = -L w /2, z c = L w /2 and 
Zd = L w /2 + hi- 

First, with the help of Eq. (17) we will derive an expression for the field in the 
far zone. The intensity distribution in this case, results to be in agreement 
with that given in [8]. Then, we will calculate the field distribution at the 
virtual source with the help of Eq. (13). Finally, Eq. (5) will allow us to find 
an expression for the field both in the near and in the far zone. We will 
the turn to analyze the case when a undulator is present, and describe the 
field from the setup as superposition of three laser-like beam from straight 
sections and undulator. In the following we will ignore the presence of the 
bending magnets, i.e. the radiation switchers. As before, we will treat them 
as if they had zero length. As a result, Eq. (95) can be simplified to 



oo) = E AB (z , r ±0 , oo) + E BC (z , r ±0 , oo) + E CD (z , r ±0 , oo) . (97) 

We already addressed the question of the applicability region of the "zero- 
length switcher approximation" in previous Section 5. 



49 



6.1 Far field from the undulator setup 



Let us describe the far field from the undulator setup in Fig. 10 by separately 
characterizing the field contributions along undulator, straight sections and 
bends and finally adding them together. 



6.2.2 Field contribution calculated along the undulator 

-> 

We first consider the contribution E BC from the undulator, following [1]. 
Assuming a planar undulator we write the following expression for the 
transverse velocity of an electron: 

cK 

v_l(z') = sin(k w z')x , (98) 



that is Eq. (52). Here K is the undulator parameter defined in Eq. (51). More- 
over, k w = 2n/A w , where A w is the undulator period, so that the undulator 
length is L = N W A W . The transverse position of the electron is therefore 

r\{z') = — cos (k w z')x . (99) 



An expression for the curvilinear abscissa s as a function of the longitudinal 
position z' is given by 

s(z') = -Lz? - sin (2k w z') , (100) 



where $ av is the time-averaged velocity along the z direction, that can be 
expressed as: 

^=4~i^)- (ioi) 



We can now substitute Eq. (99) and Eq. (98) in Eq. (17). Such substitution 
leads to a general expression, valid for any observer distance z . It is possible 
to obtain, similarly to many SR textbooks, a simplified expression valid in 
the limit for large values of z . Since we are interested in the contribution 
of the undulator device to the total field at the observer position, we will 
integrate Eq. (17) only along the undulator. Then all terms in (z - z')" 1 in 
the phase factor of Eq. (17) can be expanded around z . In the far field limit 
we can retain first order terms in z'. Dropping negligible terms we obtain 
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Ebc(z / r±.or ^) 



exp [z'O bc ] 



— sin (k w z') + 6 X 



x+O y g\, (102) 



where 



O BC = C0 



K6 X 



2c 



2c \yl x V r k ^ c 



cos(k w z') 



K 2 



sm{2k w z') 



8y 2 k w c 

The longitudinal Lorentz factor y z in Eq. (103) is defined by 



(103) 



y z = 



VI + K 2 /2 



(104) 



and is always smaller than y, because the average longitudinal velocity 
of the electron inside the undulator is smaller than that along the straight 
sections. 

In this paper we will be interested up to frequencies much lower than the 
resonance frequency i.e. co «: 2y 2 ck w 14 . Then, the contribution due to the 
term in sin^z') in Eq. (102) can always be neglected when compared with 
the maximal field magnitude of the terms in 6 X/ y. Similarly in Eq. (103), 
phase terms in cos(^„z') and sin(2fca,z') can can also be neglected. As a result, 
Eq. (102) can be simplified as 



zc 

E BC {z ,r~t ,oo) = — I dz'exp [zO BC ] [o x x + y y) (105) 
c z J 



while 



2 + eV 



(106) 



14 Although we are outside of the applicability range of the resonance approxima- 
tion, we still assume that the number of undulator period N w is large, i.e. N w » 1 
to calculate the undulator contribution. As a result, the accuracy related with the 
undulator contribution is l/(2nN w ), and not l/y 2 as in the straight section case. It 
follows that results regarding the undulator contribution to the field are not as fun- 
damental as those regarding the straight section, because the extra approximation 
N w » 1 has been used, that is specific for the undulator device. 
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Finally, note that within these approximations, the curvilinear abscissa s in 
Eq. (100) can be simplified to 



s{z') = |-z' 



(107) 



6.2.2 Field contribution calculated along the straight sections 

Let us now calculate the curvilinear abscissa along segments AB and CD. Eq. 
(107) implies that the curvilinear abscissa at the undulator edges z' = ±L w /2 
is given by ±fiL w /(2(i av ). Moreover the curvilinear abscissa must be of the 
form s(z') = z' + const., because segments AB and CD are straight lines. It 
follows that 

s ^- z, M i -j-y- z,+ w-k <io8 > 



and 



- z - — + — for z c < z < z D , 

4y 2 4y2 



(109) 



where we have used the fact that 



With the help of Eq. (17) we write the contribution from the straight line AB 
as 



: AB = ^ fdz' exp [id> AB ] (dj+ 6 y $) (111) 

C Zo J 



where O ab in Eq. (Ill) is given by 



O ab = co 



01 + e l 

2c 



-Zn ~ 



4cy 



(112) 



6 X = x/z and d y = y/z being the observation angles in the horizontal and 
vertical direction. Note that we also used y - c when calculating the term 
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L w /(4cyl) in Eq. (112). The contribution from the straight line CD is similar 
to that from the straight line AB and reads 



Ecd - 



icoe 
c 2 z 



ZD 

J dz' exp [z'Ocd] (d x x + d y yj 



zc 



(113) 



where C d in Eq. (113) is given by 

L 



2c Zo 



4cyl Acy 2 2c \y 2 



(114) 



In general, the phases C d and O^b start exhibiting oscillatory behavior 
when coz' /(2cy 2 ) ~ 1, which gives a maximal integration range in the longi- 
tudinal direction. Similarly as before, in general one has that the formation 
lengths Lf S i and Lf s2 for the straight sections AB and CD can be written as 

£-fs(i,2) ~ min[A}/ 2 ,L ( i /2) ] . (115) 
Depending on the wavelength of interest then, Lf S (i /2 ) ~ Ay 2 or Lf S (i /2 ) ~ L(i,2)- 
6.2.3 Total field and energy spectrum of radiation 

Consider Eq. (97). The contributions for segment AB and segment CD are 
given by Eq. (Ill) and Eq. (113). One obtains 



Eab — 



2e(d x x + d v y) 

(i/yi + ei + e^czo 

x < exp 



exp 



. 9 2 x + 3] 

ICO — Z 

2c 



exp 



icoL m icoL T , 



+ 



Acyl 4cy 2 



MLi+LJ2U]_ + e2 + e2 \ 



2c 



X 



exp 



.coL 



4c \y 2 

that can also be written as 

itoehi 



(116) 



Eab — 



C 2 Z 



x exp 



6>sinc 


coL\ 
Ac 


icoUv 




4c 





exp 



exp 



ito0 2 z o 
2c 



(117) 
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Similarly, 



Ecd - 



2e{d x x + d y y) 



(l/y2 + el + 0fjCZ o 

x j- exp 
+ exp 



exp 







icoL w 


icoL w 


[ l(V 2c Z °\ 


exp 


_4cy 2 z 


4cy 2 



ML 2 +LJ2)(l_ + e2 + e2 \ 



2c 



that can also be written as 



(118) 



~ icoeL 2 * . 
Eqd = — — Osmc 



C A Z n 



exp 



icod 2 z 
2c 



x exp 



icoL w 1 2 



exp 



(119) 



Finally, the contribution for the segment BC is written as in Eq. (105). Cal- 
culations yield: 



* 2e(d x x + d v y) 

Ebc = —, : exp 

(llyl + 6l + 6$cz 



. 6 X + 6 V 

ico — z 

2c 



x j - exp 
+exp 



icoL 7l 



4c 



4c 



that can be written as 



(120) 



E BC = — 5 — 0sinc 

c 2 z 



coL 



exp 



icod 2 z 
2c 



(121) 



The total field produced by the setup is obtained by summing up Eq. (116), 
Eq. (118) and Eq. (121). Note that the same definition for the observation 

angle 6 — r ± /z is used in these equations. This means that the observation 
angle is measured starting from the center of the undulator, located at z = 0. 
The energy density of radiation as a function of angles and frequencies can 
be written substituting the resultant total field in Eq. (73). We obtain 
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dcodD, n 2 c (i + yiQiy 
(l + y 2 6 2 ) 



+ exp 



. coLi 
l 2cf 



.coL 7l 



4cy 2 \ 2 



K 



2n2\ 



4cy 2 \~ 2 



1 + — + y l Q 



+ 



1/y 2 + e 2 1 
Vyl + e 2 1 

ia)L 7l 



exp 



icoL,, 



2n2\ 



4cy 2 \ 2 



1 + — +y z & 



+exp 
- exp 



K 2 



4cy 2 \ 2 
. coL 2 



1 + ^— + y 2 d 2 j 
ia)L 7 



(l + y 2 6 2 ) 



+ 



+ exp 



.coL u 



K 2 



2n2\ 



4cy 2 \ 2 



1 + — + y l Q 

1 2 



2n2\ 



'2cy 2 v 1 I 4cy 2 V ' 2 



1 + — +/0 



(122) 



that is equivalent to the analogous expression in [8]. 



6.2 Virtual source characterization and field propagation 



Expressions in Eq. (117), Eq. (119) and Eq. (121) can be interpreted as far 
field radiation from separate virtual sources. Note that the same definition 
for the observation angle 6 = r ± /z is used in these equations. Let us find the 
locations of the virtual sources. We will see that these locations corresponds 
to points along the longitudinal axis where wavefronts are plane, i.e. where 
maximal simplification arises. 

When F(0) is real, Eq. (12) describes the far field in terms of a spherical wave 
centered at z = 0. In the more generic case of a spherical wave centered 
at some position z s Eq. (12) should be substituted by Eq. (11) with z = 
z s , as the typical phase phase factor of a spherical wave centered in z s 
is exp[z'cu0 2 (z o + z s )/(2c)]. In this case F(z s ) in Eq. (11) is a real function. 
Accordingly, Eq. (13) should be modified to Eq. (14). Once we substitute Eq. 
(117), Eq. (119) or Eq. (121) in Eq. (14), there are particular values of z s (z s i 
for the segment AB, z s2 for BC and z s3 for CD) such that the phase in Q 2 is 
cancelled and that E(z s , r ± ) gives a plane wave. These are: 




(123) 



z s2 = (124) 



and 



L w L2 
z S 3 - ~2 + "2" • 



(125) 
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Substituting Eq. (117), Eq. (119) and Eq. (121) into Eq. (14) we obtain the 
following plane wavefronts describing the three virtual sources at positions 

z s i/ z S 2 and z s y. 



Eab (z s i/ r±) - — 



ico lL w L\ 




exp 



ico 

—r± 
L c 



e 



(126) 



Ebc(z S 2, r ± ) = - 



CO 



2 eL 

71 



2nc 3 



dd 6 sine 



4c 



d 2 



+ ylt 



exp 



— r, ■ 6 



(127) 



and 



~ , -> , co 2 eL 2 
E CD (z s3 ,r ± ) = -—- T exp 



x 



2nc 3 

d0 6 sine 




exp 



—r ± ■ e 

c 



(128) 



Note that L t , L 2 and L w can assume different values, y and y z are also dif- 
ferent. It may therefore seem convenient to introduce different normalized 
quantities, referring to the undulator and the straight lines. However, in 
the end we are in summing up contribution all contributions from differ- 
ent sources, so that it is important to keep a common definition of vertical 
displacement (or observation angle). Therefore we prescribe the same nor- 
malization for all quantities: 



6 = 



IcoL 



tot 



e, 



4> = 



y 2 c 



(129) 
(130) 



and 



r± = 



I co 

^totC 



(131) 



Then, we introduce parameters L\ = Li/L tot/ L 2 = L 2 /L tot , L w = L w /L tot and 
cp w = y 2 ly 2 z cp. Finally, we define z s = z s /L tot . Eq. (126), Eq. (127) and Eq. (128) 
can then be written as 
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£ab (zi, f±) =-h exp (tw^a, + Li<J») 



x (i0 sine 



exp 



if ± ■ e 



(132) 



£Bc(z2,r_ L ) = -L w I rf0 0sinc 



t(^ 2 + ^) 



exp 



/fx • 



and 



Ecd (z 3 , f±) = -U exp | - (Lw^n, + E 2 $) 

t(* + *) 



x dd 6 sine 



exp 



zfx • 9 



where we defined 



(133) 



(134) 



(135) 



Besides factors in the sine functions, the integrals in Eq. (132), Eq. (133) and 
Eq. (134) are mathematically identical to Eq. (79). One may check that, in 
the limit for <<c 1 and § w <sc 1 one obtains the same results as for edge 
radiation from a single straight section. 

Besides this obvious limit, there is a second region of interest in the param- 
eter space that can be dealt with analytically and corresponds to § <sc 1 and 
§ w » 1. In this particular limit, the contribution from the undulator, Eq. 
(133) can be neglected, because the sinc(-) is strongly suppressed. In this 
case one has the following virtual sources: 



E(zi, r ± ) = — exp 
Ei 



smc 



(l^f/Ii), 



(136) 



E(z 2 , fx) - , 



(137) 



E(z 3 , r ± ) = -t— exp 
E 2 



smc 



(|fxf/E 2 ), 



(138) 
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Note that the previous three equations describe virtual sources characterized 
by plane wavefronts. Application of the Fresnel propagation formula, Eq. 
(5) allows one to reconstruct the field both in the near and in the far region. 
Results can directly be obtained from Eq. (83) accounting for the fact that 
we want to calculate contributions at the same observation point, which 
requires substitution of z /L^ iW/2 ) with (z - z s (i,2,3))/L(i rWi 2). From Eq. (83) we 
obtain the following results for the two surviving contributions: 



E AB {z ,r 1 _) = -— exp 







4 


exp 



2(z - z sl ) 



x 



exp 
exp 



iU?\ 



\ 2(z -z sl )(L 1 +2z -2z sl ), 
2(z -z sl )(-L 1 +2z -2z sl )j 



(139) 



* -> 4f ± 
EcD(z ,r±) = -^ r exp 



iE w <p u 



exp 



2(z - z s3 ) 



x 



exp 
exp 



iL 2 f\ 



\ 2(z -z s3 )(L 2 + 2z -2z s3 ), 
2(z -z s3 )(-L 2 + 2z -2z s3 )j 



(140) 



Eq. (139) and Eq. (140) solve the field propagation problem for both the near 
and the far field in the limit for <ft <?c 1 and (p w » 1. The directivity diagram is 
obtained by summing up Eq. (139) and Eq. (140) and taking square modulus 
of the sum. 

In the particular case I4 = L 2 = L w = L tot /3 we obtain simplified expressions: 



EAB(z ,r±) - ^r- 



4f 
a~ ex P 



~2 

n 



exp 



'2(^ + 1/3) 
ir z ± 



exp 



6(z + l/3)(l+2z ), 



6(2 + l/3)(22 + l/3). 







exp 


~l2 



(141) 



E C d(z , r ± ) = -—exp 
r\ 



2(z -l/3) 



ir , 



exp I 



6(z -l/3)(2z -l/3), 
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The corresponding directivity diagram is given by 
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- exp 
+ exp 

- exp 
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exp 
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exp 
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iz 2 o 2 



6(z -l/3)(2z -l/3) / 



/z^ 



6(z -l/3)(2z -l), 



exp 



(143) 



Even though it refers to a particular case where L\ = L 2 = L w = L tot /3, 
Eq. (143) still depends on the parameter <ft w » 1. Its dependence on (p w is 
periodic, with period 127T. For the sake of exemplification, in Fig. 11 we plot 
the intensity profile for fy w = 127Z at different distances z and we compare 
these profiles with the far field asymptotic behavior. In Fig. 12 and Fig. 13 
we plot, instead, the intensity profile for different values of (p w at z = 0.6 
and in the asymptotic case for z » 1 respectively. 



7 A critical re-examination of conventional edge radiation theory 



In this Section we compare our findings with treatments of edge radiation 
and TUR, that can be found in literature. Conventional understanding is 
based on the zero-length switcher approximation. This allows direct com- 
parison of results with the contributions calculated in the previous Section 
5 and Section 6. Since there is a lot of literature dealing with edge radiation 
and TUR we will limit ourselves to a few significative works only. In par- 
ticular, we will refer to [8], [12], [9] and [7], even though other works report 
the same results. For example, the same understanding can be found in the 
very recent 15 review [13]. 

We found agreement with literature as concerns calculations in the far-field 
zone. 

As we have discussed before, the far-field asymptotic is valid at observation 
15 The year of writing is 2006. 
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Fig. 11. Evolution of the intensity profile for TUR for cp w = \2n . These profiles are 
shown as a function of angles at different observation distances z - 0.6, z - 1.0, 
z - 2.0 and z - 10.0 (solid lines). The dashed line always refers to the far-zone 
intensity. 

positions z » L, independently of (p = L/(y 2 A). It should therefore be 
stressed that in the case (p » 1 one may be in the near zone even though the 
observer is located many formation lengths y 2 A from the downstream edge 
(for example, when z = L). In the far-zone case Eq. (69) gives the electric 
field in the far region. The correspondent directivity diagram, Eq. (78), is 
plotted in Fig. 4, for different values of the main parameter of the theory, 
<ft. This result is well-known and in perfect agreement with literature. For 
example one may see that, aside for proportionality factors due to different 
choice of units, there is agreement between the radiation energy density 
in Eq. (75) and the photon flux per unit solid angle in Eq. (18) of [8] (also 
compare, qualitatively, Fig. 4 of that reference and our Fig. 4). Similarly, the 
field in Eq. (75) is in agreement with Eq. (17) of [12] or Eq. (31) of [9] (in the 
limit for a large observation distance) accounting for the fact that the origin 
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Fig. 12. Intensity profiles for TUR at zq - 0.6. at different values <p w - (l/4)(127z), 

$ w = (1/2)(127T), $ w = (3/4)(127T), ^ = (127T). 

of the reference system is chosen at the exit of the straight section by the 
author of [12] (also, the electron charge e is negative in that reference). Note 
that in the far zone, the field exhibits a spherical wavefront (independently 
of the choice of coordinate system). 

Differences between our results and results in literature arise, however, in 
the near zone. 



7.1 Near-field edge radiation 



Case for » 1. Let us pose d = z — L/2. In the case when < d «; y 2 A, we 
have seen that the intensity profile for edge radiation reproduces better and 
better (as d decreases) the single edge source, i.e. Eq. (94). According to our 
understanding then, it follows that the photon flux depends on y, does not 
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Fig. 13. Intensity profiles for TUR in the asymptotic limit for zq » 1. at different 
values 4> w = (1/4)(12ti), # w = (1/2)(12ti), # w = (3/4)(12tz), # w = (12n) 

exhibit oscillatory dependence on the transverse position r, it has a charac- 
teristic scale yA, and depends on the frequency as co 2 (Aoj /a))K 2 [cor ± /(cy)]. In 
contrast to this, the flux according to Eq. (14) in [8] 16 reads, in our notation 
(number of photons per unit surface per unit time within the bandwidth 

Aco/co) : 



dF Aoo lb 
- a — 



dS co en 2 r 2 



sm 



cor^_ 
4cd 



(144) 



where a = e 2 /(hc) is the fine structure constant and l\> is the beam current. 
This suggests that for d «: y 2 A the photon flux does not depend on y, exhibits 
oscillatory dependence on the transverse position r±_ on the scale VXd, and 
depends on the frequency as (Aa>/a>)sin 2 [cur 2 /(4ai)]. The field from which 



16 This equation coincides with Eq. (20) in [7] and Eq. (26) in the review [13]. It is a 
currently accepted expression describing the photon flux in the case under study. 
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Eq. (144) is derived is given explicitly in several papers (e.g. Eq. (24) in [9], 
Eq. (4) in [11]). In our units, it reads: 



~ -> v 4ze r ± 

E(z , r ± ) = T exp 

c r , 



.cor i 
1 4crf 



sin 



cor^ 
Acd 



(145) 



We disagree with the conclusions in Eq. (144) and Eq. (145). In particular 
the expression for the field in Eq. (145) is in contrast with our Eq. (91). 

The region of parameters for d ~ y 2 A is not discussed in literature. As we 
have seen, this corresponds to a region where the near-field contribution 
from the downstream edge is present. 

When y 2 A «: d «: L the contribution of the downstream edge is dominant, 
and the intensity distribution tends to approximate const x d 2 /(d 2 +cp) 2 , that is 
the single edge far field limit. This case is practically described in literature 
when discussing "an electron exiting a bending magnet at the upstream 
end of a straight section" of infinite length (cited from [8]. Compare, for 
example with Eq. (8) in the same reference). Although this situation is per se 
not physical, in the sense that an infinite straight section cannot be built, it 
has a practical realization (for q> » 1) in the before-mentioned observation 
region y 2 A «: d «: L. 

The region of parameters for d ~ L is not discussed in literature, because 
the near-field region is understood for d < y 2 A only. However, as we have 
seen, the case d ~ L corresponds to a region where interference begins to 
become important. In our understanding, this is a transition region between 
the near and the far zone (d » L). 



Case for <^«1. In this situation we can compare results in literature with 
our Eq. (83). In particular, Eq. (31) of [9] (or, equivalently, Eq. (17) of [12] for 

the Fourier transform of the field E) is said to be applicable at for observation 
regions where d <c y 2 A (using our notations) and for L <^c y 2 A, i.e. when 
cp <^c 1. Accounting for the fact that the origin of the reference system in [9] 
is at the downstream edge of the straight section, we can write Eq. (31) of 
[9] in our dimensional units as: 



ft ^ 2e ^ 
E(z , r±) = t" exp 

c r , 



,cor\ 




. 2CZ 0. 


exp 



icor. 



x 



expj 



2cz (l + 2z /L) 




icor. 



2cz (-l + 2z /L). 



(146) 
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Thus, Eq. (31) in [9] appears to differ of the phase factor expjz'cur 2 /[2c(L-2z )]} 
with respect to our Eq. (83), directly calculated by propagating the field at 
the virtual source position. We disagree with the result in Eq. (146) (i.e. Eq. 
(31) in [9]). The intensity pattern, however, is correct in both cases 17 . 

It is of some interest to understand the origin of the difference between our 
Eq. (83) and Eq. (31) in [9]. In this regard, in [9] it is reported that when 
L » d (and L «: y 2 A) , the interference between upstream and downstream 
edge is negligible. In this limit, Eq. (31) in [9] gives back Eq. (24) in [9] 
that is Eq. (145). On the one hand Eq. (24) in [9] is said to be valid for 
an infinitely long straight section followed by a downstream edge. On the 
other hand, Eq. (31) in [9] (and thus its limiting case for L » d as well) is 
reported to be valid for d,L <c y 2 A. Since the limiting case of Eq. (31) in [9] 
for L » d coincides with Eq. (24) in [9] one concludes that an infinitely long 
straight section is identified de facto with the case L » d. We hold this to be a 
misconception, in disagreement with our understanding. The length of the 
straight section is naturally measured by the parameter <ft = L/ (y 2 A). Such 
length is independent of the position of the observer. Only when (p » 1 it 
is possible to have radiation from a single edge. In our understanding thus, 
the condition L » d is still necessary, but not sufficient to have single-edge 
radiation. Describing the situation for (j) «: 1 we have seen from Eq. (83) 
that the total field is always given as the result of the interference of the 
two virtual sources located at the straight section edges and one never has 
radiation from single edge. In contrast to this radiation from a single edge 
(in Eq. (24) of [9]) is explicitly reported to hold for L <<c y 2 A in the limit for 
L » d. 

We conclude that the origin of the difference between our results and those 
in literature may be ascribed to a different understanding of the parameters 
of the theory. While we hold ^ to be the main parameter of our treatment, 
result in literature seem to neglect its existence. By doing this, conventional 
treatments neglect the fact that d and L should always be compared with 
the third main length scale of the problem, that is y 2 A. 



7.2 Transition Undulator Radiation 



As we have seen, TUR can be addressed as a more complicated edge- 
radiation setup. In this case we have contributions from three parts, two 
straight lines and the undulator. As we have seen, the undulator contri- 
bution is similar to a straight line contribution, the only difference being a 



Note that the intensity profile in Eq. (27) of the review [13] is the square modulus 
of Eq. (31) in [8]. 
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different average longitudinal velocity of the electron. Then, the far-zone 
region can be identified by distances z » L tot . In the far zone, well-accepted 
expressions for the TUR emission are reported in literature [14, 15, 16], that 
are equivalent to the following equation for the radiation energy density as 
a function of angle and frequency: 



dW e 2 
dcodQ n 2 c 



2y 2 dK 2 



(l+K 2 /2 + y 2 d 2 )(l + y 2 d 2 ) 



sin 



(147) 



We disagree with this expression. In our understanding, there cannot be any 
range of parameters in the setup in Fig. 10 where Eq. (147) is valid. 

In order to prove this it is sufficient to compare Eq. (147) with Eq. (122), that 
is equivalent to Eq. (20) of reference [8] 18 . In particular, since Eq. (147) does 
not depend on the straight section lengths L\ or L 2 , we conclude that the 
only regions where such comparison makes sense can be for L\ = L 2 = 0, or 

Li = L 2 — L v 



19 



When Lj = L 2 = Eq. (122) reduces to: 



dW 
dcodQ. 



n 2 c 



2y 2 6 



1 2 



sin 



1 + K 2 /2 + y 2 6 2 
that is obviously different from Eq. (147.) 
When I4 = L 2 = L w Eq. (122) becomes: 



(148) 



dW 
dcodQ. 



n 2 c 



x 



2y 2 6K 2 

(1 + K 2 /2 + y 2 6 2 ) (1 + y 2 6 2 ) 



1 




- sin 




2 


2y 2 A \ 



(1 + K'/2 + y'9*) 
¥? Sm 



1 2 



(149) 



18 We believe that there is a small misprint in Eq. (8) of [8]. The symbol "6" after 
the sign "«" should be replaced by the fine-structure constant "a" . 

19 Referring to Fig. 10, note that when L 2 » y 2 K the contribution from the straight 
section edge in D is always larger than those from the undulator edges in B and C 
at any distance in the near and in the far zone. This is the case because outside the 
undulator the longitudinal velocity is nearer to c than inside {y\ < y 2 ). Thus, for 
long straight sections L\ = L 2/ with L 2 » y 2 A, the straight section contribution will 
be dominant compared to the undulator contribution. 
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that is also different from Eq. (147) 20 . We conclude that Eq. (147) has no 
physical meaning. 

Finally, we put attention on the fact that there is a general need, in the FEL 
community, to extend the current theory of TUR to cover the near zone. 
Some attempt in this direction has been reported in [20]. The author of [20] 
discusses a possible use of coherent TUR to produce visible-to-infrared light 
synchronized with X-rays from an X-ray free-electron laser. For simplicity 
we will ignore the delicate issue of technical realization of the scheme pro- 
posed in [20] for the LCLS facility, and we will restrict our discussion of a 
issue pertaining fundamental questions of electrodynamics. 

In our view, the extension of the theory of TUR in Appendix A of [20] to the 
near zone also includes some misconceptions. The treatment begins with Eq. 
(Al), that is an expression for the vector potential A at a given observation 
position and at a given frequency of interest. We report this expression here, 
for the readers' commodity: 



A(co) 



T/2 

I 

-T/2 



exp[ico(t + R(t)/c)]- 



ft x (ft x /?) 



(150) 



Here T = L/c is the time that a photon takes to travel the undulator length, 
R(t) is the distance from the electron position at retarded time t to the 
observer, the unit vector ft is the direction from the retarded position of 

the electron to the observer and jS(f) is the electron velocity in units of the 
speed of light. Eq. (150) is a modification of the well-known Eq. (14.62) of 
[21]. In fact, integration limits in Eq. (150) have been changed from ±oo to 
+T/2 "under the assumption that the electrons go on straight paths before 
and after the undulator" (cited from [20]). As a first remark, as discussed 
above, the previous assumption is not justified in our view, as it does not 
make sense to discuss generically about infinite straight lines preceding and 
following the undulator. However, we believe that there is another problem 
with the use made of Eq. (150). The author of [20] starts with Eq. (150) and 
modifies it accounting for the fact that, in the near zone, ft is not constant. 
The problem with this is that the derivation of Eq. (150) relies a priori on the 
assumption that ft = const. In particular, the classical result 



rtx[(rt-p)xp\ ld_ 
{l-n-$)\r -r'{t)\ r odt 



ft x (ft x p) 



(151) 



20 For example, in the limit for K «: 1 Eq. (147) tends to zero, while Eq. (149) gives 
back Eq. (75) with L = 3L W as it should be. 
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is used, where r = \r \ is the distance of the observer from the origin of the 
coordinate system. Eq. (151) assumes n = const. 

Our conclusion is that Eq. (150) cannot be manipulated assuming that n is 
not constant, as has been done in [20]. When n is not constant, as in this 
case, Eq. (151) should be modified to the expression presented before in Eq. 
(46) that may be found in the interesting but perhaps little-known reference 
[23], dating back more than twenty years. 

However, in this paper we have seen that it is possible to develop a more 
convenient theory of near-field Synchrotron Radiation Theory (including 
the case of TUR) in the space-frequency domain, where Fourier Optics 
techniques can be taken advantage of. 



8 Conclusions 

In this paper we presented a general connection between Fourier Optics 
and classical relativistic electrodynamics. Consistent use of such connection 
resulted in the formulation of a modern theory of Synchrotron Radiation in 
terms of laser-beam optics. 

In particular, we developed a theory of near-field SR in the space-frequency 
domain based of Fourier Optics techniques. These techniques can be taken 
advantage of without limitations for SR setups, because the paraxial approx- 
imation can always be applied in the case of electrons in ultra-relativistic 
motion. We restricted ourselves to the analysis of single-particle radiation. 
As we demonstrated in Section 2, and discussed in Section 3 the use of 
paraxial approximation allows reconstruction of the field in the near-zone 
from the knowledge of far-field data only. The solvability of the inverse 
problem for the field starting from far-field data allows characterization of 
any synchrotron radiation setup by means of virtual sources. In cases of in- 
terest these sources exhibit a plane wavefront, and can be pictured as waists 
of laser-like beams. These laser-like beams help developing our theory in 
close relation with laser-beam optics. In particular, usual Fourier Optics 
can be used to describe the field at any distance, thus providing a tool for 
designing and analyzing SR setups. 

We gave several examples of applications of our theory in Section 4, Section 
5 and Section 6. These examples provide both physical insight and under- 
standing of several non-intuitive phenomena. In Section 4 we treated the 
case of undulator radiation. In the case of perfect resonance we found an 
analytic expression for the virtual source and we propagated such expres- 
sion both in the near and in the far zone. The applications in the following 
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Section 5 and Section 6 actually constitute the first comprehensive theory of 
edge radiation and transition undulator radiation. 

In Section 5 we dealt with the fundamental case of an electron moving on 
a straight section. Such kind of system can be considered, in all generality 
as a fundamental building block for a number of more complicated setups. 
We demonstrated that two equivalent pictures in terms of laser-like sources 
can be presented. The first takes advantage of a single virtual source in the 
middle of the straight section. The second relies on two virtual sources at the 
edges. The virtual sources at the edges can always be described analytically 
independently of the region of parameters considered. The main parameter 
of the theory is identified to be the ratio between the straight section length 
L and y 2 A. When this ratio is small the single-source picture is found to be 
natural, and two asymptotic regions of observation are found, the near and 
the far zone. When the ratio is large, the two-source picture is more natural 
and one is left with four asymptotic regions of observation: far zone for 
two interfering sources, near zone for two interfering sources, far zone for 
a single edge and near zone for a single edge. Results in Section 5 can be 
directly used to describe an edge-radiation setup in the case upstream and 
downstream bends can be considered to have zero length. These bends act 
like switchers, switching the radiation harmonic (and the beam current har- 
monic) on and off. Thus, our results for the straight section can be applied to 
any situation with zero-length switchers. Our expressions are of importance 
also in the case when switchers cannot be considered having zero-length. 
In fact in general, as said before, straight sections can be considered as basic 
building blocks for any magnetic setup. 

To demonstrate this last point, in Section 6 we used results of Section 5 as 
building blocks for describing the case of TUR in the zero-length switcher 
approximation. In particular, we dealt with a setup composed by an un- 
dulator preceded and followed by two straight sections. In the zero-length 
switchers approximation, such setup describes TUR emission. Although 
more sophisticated with respect to the straight-section motion, the example 
in Section 6 does not present any conceptual difficulty in addition to those 
encountered in the treatment of the straight-section motion. 

Finally, in Section 7, we compared our findings with current understanding 
of edge radiation and TUR. We found that several misconceptions exist in 
the conventional understanding of edge radiation in the near zone, as well as 
in the treatment of TUR. We extensively discussed these misconceptions and 
proposed our understanding of the correct description of these phenomena. 

We find it interesting to remark that, although classical relativistic electro- 
dynamics is a relatively old branch of physics, there are many novel and 
valuable results pertaining this field that have been obtained only recently. 



68 



Undoubtedly, many more will follow. These recent developments illustrate 
a recurring fact in the history of Physics. Namely, subjects appearing well- 
understood, and perhaps even a little old-fashioned may still have surprises 
in store for us. 
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